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PREFACE 


This  document.  No.  FAA-rd-^O-JC;  is  the  final  report  of  a  study  conducted 
by  the  Gates  Learjet  Corporation  to  develop  an  improved  and  simplified 
Procedure  for  predicting  and  evaluating  engine  inlet  anti-icing  system. 

All  work  was  performed  in  the  Advanced  Design,  Aerodynamics,  and  Propuls¬ 
ion  Analysis  departments,  under  the  coordination  of  Mr.  A.  M.  Heinrich, 
Program  Manager.  The  technical  analysis  was  performed  by  Messrs.  R.  Ross, 
N.  Ganesan,  D.W.  Newtor  and  R.  Sundccfst.  Development  of  the  anti¬ 
icing  analysis  computer  programs  v  <s  originally  formulated  by  Mr.  T.  M. 
Kutty  prior  to  his  departure  from  the  company 

Appreciation  is  fatefully  extended  for  the  cooperation  and  assistance 
provided  bv  r.rsonnel  of  the  NASA-Lewis  Research  Center  in  the  conduct  of 
icing  tunnel  tests. 
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l'.O  INTRODUCTION 


Gates  Learjet  introduced  two  new  aircraft,  Models  35  and  36,  into  their 
product  line  during  1974.  These  aircraft  are  powered  by  Garrett-Ai Research 
TFE  731-2  turbofan  engines  and  use  bleed-air  for  thermally  anti-icing 
engine  nacelle  inlet  lips,  wings,  and  horizontal  tail.  During  development 
of  these  aircraft  it  was  founu  that  existing  techniques  for  predicting  anti¬ 
icing  system  performance  were  either  very  complex  or,  when  simplifying 
assumptions  were  made,  were  highly  questionable  in  accuracy.  In  addition, 
current  procedures  for  anti-icing  system  certification  include  a  significant 
amount  of  flight  testing  in  simulated  and  difficult  to  find  natural  icing 
conditions  which  is  an  expensive  and  time  consuming  process  providing  a 
very  limited  sample  of  the  actual  icing  environment.  Also,  the  validity 
of  extrapolating  data  obtained  from  such  flight  tests  to  other  FAR  25 
Appendix  C  conditions  has  always  remained  questionable. 

What  has  appeared  to  be  needed  is  a  simpler,  less  costly,  and  less  time 
consuming  technique  where  flight  test  requirements  can  be  reduced  to  a 
minimum.  To  achieve  this  a  procedure  is  needed  which  would  utilize 
analytical  prediction  methods  to  evaluate  the  full  scope  of  the  icing 
environment  and  icing  wind  tunnel  tests  to  confirm  the  method.  The 
procedure  must  be  sufficiently  accurate  and  comprehensive  to  be  accep- 
able  to  the  FAA,  but  at  the  same  time  be  simple  and  economical  enough  to  be 
widely  acceptable  within  the  industry.  J _ 

Against  this  background  but  focusing  on  only  one  area  of  a  transport  type 
airplane,  a  contract.  No.  DOT  FA76WA-3852,  was  granted  by  the  FAA,, 
Washington,  D.C.,  to  develop  a  procedure  for  predicting  and  evaluating  the 
performance  of  engine  inlet  anti -icing  systems  for  compliance  ,.ith  FAR  25 
ice  protection  requirements.  The  desire  was  to  minimize  the  need  for 
conducting  flight  tests  in  natural  icing  conditions.  The  procedure  was  to 
include  consideration  of  water  droplet  impingment  and  collection 
efficiency,  internal  and  external  heat  transfer^  and  mass  transfer  of  the 
impinging  water.  A  large  spectrum  of  environmental  and  operational 
factors  were  to  be  analyzed  in  development  of  the  procedure  to  identify 
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those  environmer tal/operational  conditions  which  should  be  recommended  for 
evaluation  by  the  procedure.  Icing  wind  tunnel  tests  were  required  to 
assist  in  the  procedure  development.  ■ 

The  following  report  describes  the  method  of  approach  taken,  icing  wind 
tunnel  tests  conducted,  and  a  correlation  between  predicted  results,  using 
the  methodology,  and  the  tunnel  test  results.  Conclusions  and  recommen¬ 
dations  are  presented  plus  details  of  the  computer  programs  developed. 


2.0  ANALYSIS  METHODOLOGY 

Prediction  of  anti -icing  system  performance  is  based  on  a  combination  of 
aero-thermodynamic  theory  and  empirical  relationships.  A  summary  of  these 
methods,  based  on  the  works  of  many  investigators,  was  prepared  and  pub¬ 
lished  by  the  Federal  Aviation  Administration  in  1964)  Solution  of  the 
general  problem  of  internal  heat  transfer  and  of  heat  and  mass  transfer 
from  a  wetted  surface  in  forced  convection  is  quite  involved  and  tedious. 
The  method  of  solution  customarily  involves  several  trial -and-error 
calculations  that  are  intermediate  between  the  final  answer  and  the  basic 
factors  that  define  a  particular  anti-icing  situation.  In  order  to 
lessen  the  burden  of  this  effort  a  digital  computer  program  was  developed 
to  perform  these  calculations. 

The  solution  desired  in  anti-icing  calculations  is  for  the  internal  heat 
transfer  rate  of  thermal  heat  systems,  the  external  heat  transfer  by 
convection,  evaporation  of  surface  water,  and  the  sensible  heat  change  of 
the  impinging  water.  After  determining  the  heat  transfer  characteristics, 
the  anti-icing  performance  (i.e.,  the  rate  of  water  impinging  on  the 
surface  and  the  rate  of  evaporation)  can  be  calculated  for  any  known  set 
of  flight  and  atmospheric  conditions.  The  difference  between  the  rate  of 
water  impinging  on  the  surface  and  the  rate  it  is  evaporated  is  the  run- 
back  rate.  Knowing  the  run-back  rate,  the  amount  of  run-back  ice  corres- 
pounding  to  the  time  taken  to  travel  through  the  horizontal  extent  of  the 
cloud  is  calculated. 

The  analysis  is  made  using  a  poi.nt-by-point  approach  by  dividing  the 
heated  area  into  several  small  segments.  This  method  represents  a  more 
realistic  evaluation  because  skin  temperature  and  impingement  rate 
gradients  are  generally  too  steep  for  use  of  average  values. 


1.  Bowden,  D.T.,  et.  al.:  "Engineering  Summary  of  Airframe  Icing  Technical 
Data".  Technical  Report  ADS-4,  Federal  Aviation  Administration,  March 
1964. 
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2.1  Aero-Thermodynamic  Relationships 

The  theoretical  and  empirical  relationships  used  in  calculating  anti¬ 
icing  performance  of  a  system  are  arranged  into  several  categories  for 
convenience.  The  grouping  is  in  a  sequential  order  of  calculation; 

i.e.,  the  equations  presented  in  a  particular  group  use  the  known 
parameters  mentioned  in  that  group  plus  those  that  are  mentioned  (known 
or  calculated)  in  preceding  groups.  Various  parameters  are  identified 
by  symbols  which  are  identical  or  similar  to  the  variables  used  in  the 
computer  program.  This  helps  to  easily  understand  the  logical  development 
of  the  computer  programs  from  a  listing  of  their  statements. 


Basic-Relationships 


Known  parameters: 

Altitude  (ALT),  ft 

Ambient  Temperature  (T^),  °R 

True  Airspeed  (VQ),  ft/sec 

Gas  Constant,  R  *  53.35  (ft-lb)/(lb- °R) 

Based  on  the  above,  the  following  equations  can  be  developed 

Specific  heat2,  C  =  .2355  +  7.6  x  10"6  T.  Btu/(lb-°R)  (2-1) 
P  " 

Coefficient  of  viscosity* 


7.475  x  ICfV-5 

MU  *  - - - -  lb/(ft-sec)  (2-2) 

TA  +  216 

Coefficient  of  Thermal  Conductivity 3 

Kj  =  (.06944  Ta  +  4.722)  x  10-7  Btu/(ft-sec-°R)  (2-3) 


2.  Eshbach,  O.W.  and  Souders,  M.:  "Handbook  of  Engineering  Fundamentals". 
Third  Edition,  Wiley  and  Sons,  New  York,  1974,  (P.843,  Table  3). 

2.  ibid,  (P.576,  Figure  4) 

3.  Baumeister,  T.:  "Standard  Handbook  for  Mechanical  Engineers". 

7th  Edition,  McGraw-Hill,  New  York,  1966,  (P.4-93,  Table  2). 
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Prandtl  number,  Pr  =  ( Cp )  (MU)/K^ 

Ambient  pressure 

P  =  2116.21/  e[ALT/(2771°  '  *098774  ALT):1 
Density,  RHO  =  PA/(R  T^) 

Reynolds  Number  per  foot,  Rg  =  (RHO)  (Vq)/(MU) 

The  equations  for  specific  heat,  coefficient  of  viscosity  and  coefficient 
of  thermal  conductivity  were  developed  by  a  curve  fit  to  the  data  presented 
in  the  cited  references. 

the  ICEOFF  computer  routine  was  written  to  calculate  some  basic  aircraft 
performance  parameters  and,  by  calling  an  engine  performance  subroutine, 
determined  extracted  bleed  air  data  for  the  given  conditions.  Since 
engine  performance  subroutines  may  not  always  be  available  and  the  desired 
data  is  otherwise  known,  the  engine  subroutine  may  be  bypassed  and  bleed 
air  data  read  in  directly. 

Thrust/Bleed-Air  Calculations 
Known  parameters: 

Airplane  weight  (W),  lb 

2 

Airplane  wing  reference  area  (Sw),  ft 
Profile  drag  coefficient  (Cqq) 

Induced  drag  coefficient  factor  (DCD) 

Nacelle  highlight  diameter  (DIA),  ft 
Based  on  the  above,  the  following  equations  can  be  written: 


Lift  coefficient,  CL  =  W/C^j|  V2  SyD 

(2-8) 

2 

Drag  coefficient,  C.  s  C„  +  (DCD)  C, 

.  uo  L 

(2-9) 

d«9.  0  -  CD  5$  V*  Sw  lbs 

(2-10) 

Thrust  per  engine,  FN  a  D/2  lbs  (for  2  engines) 

(2-11) 

Knowing  ALT,  VQ  and  Tft,  the  bleed-air  temperature  (Tg)  and  mass  flow 
rate  (Wg)  can  be  determined  from  engine  performance  data. 


(2-4) 

lb/ft2  (2-5) 

1 b/f t3  (2-6) 

1/ft  (2-7) 


Impingement  Calculations 

The  inertia  parameter,  Kg,  the  Droplet  Reynolds  Number,  R^,  and  the 
water  collection  efficiency,  E^,  are  required  to  completely  determine 
the  rate  of  water  catch  and  the  region  of  impingement  on  a  given  airfoil 
shape. 

Known  parameters: 

Chord  (nacelle  length)  (C),  ft 
Thickness  (t),  ft 

2 

Heated  surface  area  per  foot  of  circumference  (S^),  ft  /ft 
Droplet  diameter  (D^),  ft 
Liquid  water  content  (LWC),  g/nf* 

2 

Acceleration  due  to  gravity,  6  *  32.174  ft/sec 
Surface  Pressure  coefficient  (Pc) 

Figure  2-1  illustrates  this  geometry  definition.  Based  on  the  above,  the 

following  relationships  can  be  written: 

Droplet  Reynolds  number,  R  =  (R  )Dn  (2-12) 

eD  eft  u 

Droplet  range  ratio  (LOLS):  droplet  range  is  the  distance  the 
drop  of  water  would  travel  before  impingement  if  projected  into 
still  air  with  a  given  velocity 


L0,5  n  Droplet  range  projected  into  still  air _ 

Droplet  range  projected  into  still  air  per  Stoke 's  Law 

Physically  this  is  on  average  value  of  drag  force  times  Rffl  for  a  drop 

eD 

projected  into  still  air. 


LOLS  *  .93  -  .134  In  (R  )1 


Reo  *200 
R_  >200 


LOLS  »  .74  -  .0887  In  (Ra  )  RA 

eD  B0 

These  relations  for  LOLS  were  obtained  by  a  curve  fit  of 

information  presented  by  Bowden1. 

Inertia  Parameter,  1^ 


.108  ——  0ZQ 


1.  op.  cit. 


C(MU) 


(2-13) 

(2-14) 


(2-15) 
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FIGURE  2rl  Nacelle  Configuration  Definition 
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K0  is  obtained  by  a  curve  fit  of  information  presented  by 

i 

Bowden  ' 

Impingement  parameter,  KQ  =  1^  ( LOLS )  (2-16) 

Kq  is  a  parameter  introduced  by  Langmuir  and  Blodgett1*  so  that  the 
collection  efficiency,  EM,  versus  the  inertia  parameter  curves  for  various 
droplet  Reynold's  numbers  may  be  collapsed  into  a  single  curve  for  bodies 
of  the  same  geometrical  shape.  Collection  efficiency  is  given  by: 


em  =  0 

Ko< 

.004 

(2-17) 

Em  =  .0873  [5.522  +  In  KQ3 

v; 

o 

o 

Ko< 

.01 

(2-18) 

Em  =  .08  +  .31  [2  +  .4342  In  KQ3 

.01  s 

*0* 

.4 

(2-19) 

Em  is  the  ratio  of  the  amount  of  water  intercepted  by  the  airfoil  to  the 
amount  of  water  contained  in  the  volume  of  cloud  swept  out  by  the  airfoil. 
The  expressions  shown  for  EM  are  obtained  by  a  curve  fit  of  information 
presented  by  Bowden1. 


Water  catch. 


WM  =  .623  VQ  (LWC)  JcEj,  10"4 

lb/ sec  -  ft 

(2-20) 

Water  collection  rate  (average). 

\v  ■  VSH 

lb/(sec-fr) 

(2-21) 

Water  collection  rate  at  leading  edge. 

Mm  *  2  H. 
lE  wA  V 

lb/(sec-fc  ) 

(2-22) 

1.  op.  cit. 

4.  Lagmuir,  I.  and  Blodgett,  K.:  "A  Mathematical  Investigation  of  Water 
Droplet  Trajectories".  AAFTR  5418,  Feb.  12,, 1946 

1.  op.  cit.  (Figure  2-9) 
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Water  collection  rate  at  any  point  cine  to  impingement. 

Mu  -  Mu  C1-.385(3S,)1*7S:  0  <  S'sl/3  lb/(sec-ft2) (2-23) 

"IMP  WLE 

M,.  =  1.177  My  C(l-S')1*.6:  1/3  <  S’<1  lb/(sec-ft2)  (2-24) 

"IMP  WLE 

Where  S  =  distance  (ft)  along  surface  from  stagnation  point  to  the  point 

under  consideration,  S'  =  S/S^  and  =  S^/2  (S^  =  Distance  along  surface 

from  stagnation  point  to  limit  of  impingement,  ft).  These  expressions  for 

Mu  were  obtained  by  a  curve  fit  of  data  presented  by  Neel5. 

"IMP 

Heat  Transfer  Calculations 

Heat  transfer  calculations  are  based  on  the  fact  that  at  any  point  on 
the  surface,  the  skin  temperature  attains  a  steady  value  which  allows 
equilibrium  between  internal  and  external  heat  flows.  The  formulation  of 
the  problem  is  based  essentially  on  the  development  presented  by  Gray6 
and  Gelder7.  The  equations  are  quite  lengthy  and  a  number  of  simplifying 
assumptions  and  considerable  manipulation  reduce  them  to  a  manageable 
level  as  described  in  the  following  material. 

Heat  transfer  in  forced  convection  from  a  surface  subjected  to  water 
impingement  is  represented  by  the  total  heat  transferred  by  convection, 
evaporation,  and  the  sensible  heat  change  of  the  Impinging  water.  An 
accounting  is  made  for  heat  generated  by  friction  and  the  kinetic  energy 
of  the  water  droplets.  The  expression  for  heat  output  can  be  written 
in  several  ways  depending  on  how  the  various  terms  are  grouped.  Following 
the  development  of  Gelder7: 


5.  Neel.  C.G.:  "A  Procedure  for  the  Design  of  Air-Heated  Ice  Prevention". 
NACA  TN  3130,  1954.  (Figure  16a) 

6.  Gray,  V.H.:  "Simple  Graphical  Solution  on  Heat-Transfer  and  Evaporation 
from  Surface  Heated  to  prevent  Icing".  NACA  TN  2799,  1952. 

7.  Gelder,  F.P.,  et.  al.:  "Icing  Protection  for  a  Turbojet  Transport 

Airplane:  Heating  Requirements,  Methods  of  Protection,  and  Performance 
Penalties".  NACA  TN  2866,  1953.  " 
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In  order  to  eliminate  the  lengthy  calculation  involved  in  the  solution  of 
the  equation  (2-2b),  the  following  set  of  assumptions  were  made: 

a.  The  flow  over  the  body  being  almost  adiabatic,  the  local 
stream  vapor  pressure  can  be  represented  by6: 

eL  =  eA  PL/PA  lb/ft2  (2'30) 

b.  The  local  velocity  and  pressure  may  be  belated  by  the  incom¬ 
pressible  dry  air  relation6: 

V  *  [Vo  -  iwA  <PL  -  PA>]  «/**  ..  <M1> 

c.  The  exponent  on  the  Prandt  number  is  chosen  as  .5  corres¬ 
ponding  to  the  conservative  case  of  laminar  flow  rather  than 
1/3  corresponding  to  the  case  of  turbulent  flow. 

d.  Radiation  is  small  and  is  neglected. 

e.  Conduction  along  the  skin  is  neglected  because  the  thickness 
is  small  and  the  conduction  is  small  compared  to  convection. 

f.  The  specific  heat  of  water  at  constant  pressure,  Cp  is 

1  Btu/(1 b  -°R)  .  w 

With  the  use  of  the  above  assumptions,  the  equations  can  be  combined  and 
reduced  to  the  final  form7: 


^OUT  =  H  SH  i 

!  <Ts  •  V  »  +  V’ 

■ 

p;5  +  Mw  Cp/H  )  /2GJCp 

, 

+  0.622 

Lkw  IH  eA  \ 

it  [rL  -rj 

+  rta 

(i  -  p;5)(i-pl/paK  J 

Btu/sec-ft 

(2-32) 

6.  op.  clt* 

7.  op.  cit. 


2-9 


From  (2-22)  the  following  temperature  terms  can  be  obtained6 

T1  *  <Ts  -  V  (1  +  VH)  0R 

T2  =  Vo2  (Pr5  +  MWCP/H)/2GJ  cp  °R  . 


(2-23) 


(2-34) 


I3  =0.622  L(es/PL)/Cp 


(2-35) 


T4  =  0.622  L(eA/PA)/Cp 


(2-36) 


VRTA  d-P;5)(l-VPA^JCP 


(2-37) 


Equation  (2-32)  can  then  be  written  as: 


Q0UT  =  H 


SH  T1  '  T2  +  ^  (W  +  T5  Btu/sec-ft 

m 


(2-38) 


The  analytical  procedure  developed  by  Gray6  defines  the  terms  and 
as  follows: 


T3  =  2760  es/PL 

T4  *  2'60  VPA 


(2-39) 

(2-40) 


Where  L  and  Cp  were  assumed  to  be  1066  Btu/lb  and  0.24  Btu/(lb-°R) 
respectively  and  es>  eA  are  partial  pressures  of  water  vapor  (corresponding 
to  saturated  air  unless  otherwise  noted),  at  surface  and  ambient  conditions, 
respectively.  Curves  of  these  values,  good  for  a  wide  range  of  temperature, 
were  presented  by  Gray6.  However  for  computer  use,  a  curve  fit  to  a 
psychrometrlc  chart  vapor  pressure  line  presented  in  Eshbach2  Is  used  for 


es  and  eA. 


e^,  eA  ■  144  exp 


(A  +  BZ  +  CZ2  +  DZ3  +  EZ4)  lb/ft2 


(2-41) 


Where: 


exp  Is  the  base  of  the  Napierian  logarithm  system 


-,  where  T  Is  static  temperature,  °R 


6.  op.  cit. 
2.  op.  cit. 
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If:  T  <  491.688 
,Ice 


If;  491 .688cTs671 .688 
Water 


If:  T>671.688 
Vapor 


A  =  +19.598997 
B  =  -10.431025 
C  =  -0.27550673 
D  =  +0.039^04393 
E  =  0 


A  =  +13.435296 
B  =  -5.0988424 
C  *  -1.6896174 
D  =  +0.17829154 
E  =  0 


A  =  +16.825544 
B  =  -14.213106 
C  =  +7.5567694 
D  =  -4.0151569 
E  =  +0.71697364 


Additional  parameters  that  are  required  for  completing  the  calculations 
are  latent  heat  of  vaporization,  L,  evaporation  fraction,  K,  elemental 
runback  rate,  DRB,  water  collection  rate,  end  the  change  in  bleed  air 
temperature  along  the  heated  surface. 


The  latent  heat  of  vaporization  for  water  is  represented  by: 

L  *  1348.21  -0.5620  T$  460°R  <;  Ts  £  C60°R  Btu/lb  (2-42) 

This  relation  was  developed  by  a  curve  fit  to  information  prepared  by 
the  SAE8. 

The  evaporation  fraction  is  defined  as: 

K  a  (T3  -  T4)  H/L  Mw  (2-43) 

The  value,  of  K  depends  on  equilibrium  conditions  and  i$  solved  by  iteration 
during  the  calculation. 

The  elemental  runback  rate  is: 

DRB  =  (1-K)  Mw  lb/(sec-ft2)  (2-44) 

For  the  succeeding  element  the  water  collection  rate  is: 

-  M  (at  thit  point)  +  DRB  lb/(sec-ft2)  (2-45) 

w  WIMP 

The  drop  in  bleed  air  temperature  is: 

dtb'W(ubcp>  ,R  <2-46> 

D 

this  leads  to  the  necessity  of  calculating  the  heat  transfer  coefficients. 


8.  "SAE  Aerospace  Applied  Thermodynamics  Manual".  Second  Edition,  Society 
of  Automotive  Engineers,  Inc.,  New  York,  Oct.  1969.  (f\  165,  Fig.  2C-1). 
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Internal  Heat  Transfer 


Known  parameters: 

Bleed-air  flow  rate  (Wg),  Ib/sec-ft 

Bleed-air,  temperature  (Tg),  BR 

Skin  temperature  (Tg),  °R 

Local  channel  efficiency  (EFF) 

Based  on  the  above  the  following  .quations  can  be  written: 

Specific  heat  of  bleed-ai-2 , 

CD  =  .2365  +  7.6  x  10‘6  TB  Btu/(lb-cR) 

.  PB  B 

Heat  input,  QIN  -  Wg(Cp  )  (Tg-Tq )  Btu/sec-ft 

B  . 

External  Heat  Transfer 
Known  parameters: 

Leading-edge  radius  (R^), 

Local  pressure  coefficient  (P  ) 

w 

Mechanical  equivalent  of  heat,  J  =  778  ft-lb/Btu 
Tne  following  relationships  can  now  be  developed: 

Reynolds  number  in  the  leading -edge  region, 

ReLE  “  ReftRLE 

Local  velocity,  *  VQ  (1-PC)'^ 

Reynolds  number  at  any  point,  R^  *  Re^t  S  V^/VQ 


(2-47) 

(2-48) 


(2-49) 

(2-50) 

(2-51) 


The  dry  air  convectivq  heat  transfer  coefficient  is  determined  by  the 
method  of  Gelder7.  In  the  region  of  the  stagnation  point,  an  empirical 
equation  for  a  cylinder  Is  used: 

H  -  \o.S7Kff  ReJ/\E  j  J*  1  -  (e/90)3]  Btu/(sec-ft2-°R)  (2-52) 


2.  op.  cit. 
,7.  op.  cit. 


where  0  is  the  angle  in  degrees  from  the  stagnation  point  to  the  point  of 
interest.  This  equation  is  used  up  to  Q  =  25  degrees.  Beyond  G  =  25 
degrees,  a  flat  plate  laminar  flow  equation  is  used: 


H  *  0.332  K.  P1/3 
1  r 


Pj/3  Rej5/S 


Btu/(sec-ft  -  R)  (2-53) 


where  S  is  the  distance  along  the  surface  from  the  stagnation  point  to  the 
point  of  interest.  This  equation  is  used  up  to  a  Reynolds  number  of 
2  x  10®.  From  Re^  =  2  x  10®  up  to  the  point, where  Re<-  *  1.2  x  10®  is 
treated  as  a  transition  region  and  for  Re^  greater  than  1.2  x  10®  is 
considered  fully  turbulent  flow.  For  a  flat  plate  in  turbulent  flow: 


H  =  0.0296  K 


P_1/3  R^u/S 


.8, 

'e$ 


Btu/(sec-ft^-°R)  (2-54) 


For  the  transition  region  a  linear  variation  of  the  heat  transfer 
coefficient  with  Reynolds  number  between  laminar  flow  and  turbulent  flow 
is  used.  If  the  computed  Reynolds  number  indicates  a  transitiu.  flow 
region,  the  following  procedure  is  used. 


1. 


2. 


Determine  the  position  on  the  surface,  Sy^,  where 
Re^  =  2  x  10®  and  a  compute  a  laminar  flow  heat  transfer 
coefficient  for  conditions  at  that  point  using  equation 
(2-53),  Hlam 

Determine  the  position  on  the  surface,  SyURg  where 

Re<j  =  1.2  x  10®  and  compute  a  turbulent  flow  heat  transfer 

coefficient  for  conditions  at  that  point  using  equation  (2-54), 

hturb 

Calculate  the  heat  transfer  coefficient  for  the  point  of 
interest  from: 


H  =  H 


hturb  -  hlam 


LAM  +  r7i^ 


(R  -  2  X  105)  Btu/(sec-ft  -°R)  (2-55) 


-s 
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The  local  static  pressure  is  determined  from  the  incompressible  relation. 


P.  =  P.  +  RHO  V2  P  lb/ft2  (2-56) 

L  ”  2G  0  c 

Channel  Efficiency 

Dry  air  test  data  is  used  in  determining  an  effective  anti-icing 
efficiency.  This  is  a  measure  of  the  effectiveness  of  the  hot  bleed  air 
in  heating  the  external  surface.  For  dry  air,  Mw  and  Kw  are  zero  a(nd 
(2-38)  reduces  to: 

Qout  =  H  SH  (Tl  *  T2  +  V  Btu/sec-ft  (2-57) 

and 

Tj  =  Ts  -  Ta  °R  (2-58) 

T2  =  V2  P^5/2GJ  Cp  °R  (2-59) 

Tg  is  the  same  as  given  in  (2-37). 

The  heat  supplied,  QIN,  can  be  computed  from  (2-48)  and. the  channel 
efficiency  is  then  given  by: 


Eff  =  Q  JO. 
your  xin 


(2-60) 


2.2  COMPUTER  PROGRAMS 

The  Channel  Efficiency  Computer  program,  CHANEFF,  is  described  in  Appendix 
A.  A  flow  chart,  a  program  listing,  an  example  of  the  input  format  and 
an  example  of  the  output  are  presented.  In  a  similar  fashion  the  icing 
analysis  computer  program,  ICEOFF,  is  described  in  Appendix  B. 

Included  in  Appendix  A  is  a  listing  of  a  function  subprogram,  TRP,  that 
is  called  by  both  CHANEFF  and  ICEOFF  to  interpolate  the  input  data  arrays 
for  values  of  pressure  coefficient,  surface  temperature  and  channel 
efficiency  intermediate  to  those  supplied.  Included  in  Appendix  B,  is  a 
listing  of  a  function  subprogram,  PP,  that  is  used  to  compute  the  partial 
pressure  of  water  vapor  at  saturation  temperature. 


3.0  ICING  WIND  TUNNEL  TESTS 

Icing  wind  tunnel  tests  were  conducted  to  provide  experimental  data  that 
would  assist  in  the  development  of  evaluation  procedures  for  engine  inlet 
anti -icing  systems.  The  scope  of  testing  was  selected  to  cover  a  wide 
spectrum  of  environmental  conditions  that  were  practicable  in  the  wind 
tunnel  described  in  the  following  section. 

3.1  Mind  Tunnel  Facility 

Through  sponsorship  by  the  Federal  Aviation  Administration,  approval  was 
obtained  for  use  of  the  NASA  Lewis  Research  Center  Icing  Research  Tunnel 
(IRT)  located  at  Cleveland,  Ohio.  This  facility  is  a  closed  return  tunnel 
having  a  6  x  9  foot  (1.83  M  x  2.74  M)  test  section.  Performance  capabili¬ 
ties  are  outlined  in  the  following  table: 

Airspeed:  260  KTAS  (483  KM/Hr)  maximum  without  model  installed 
Air  Temperature:  -22°F(-30°C)  to  32°F  (0°C) 

3 

Liquid  Water  Content:  Approximately  0.5  to  2.0  g/m 
Median  Droplet  Size:  11  to  20  um 

Airspeed  capability  with  the  TFE731  nacelle  installed  was  approximately 
240  KTAS. 

Air  temperature  in  the  tunnel  was  measured  by  three  (3)  total  temperature 
thermocouples  calibrated  by  NASA-Lewis  personnel  and  located  at  the  upstream 
corner  from  the  test  section.  These  were  located  on  cross-bars  attached  to 
the  turning  vanes  and  sensed  air  temperature  at  the  top,  center,  and  bottom 
in  the  tunnel  cross  section.  Tunnel  temperatures  were  monitored  by  both 
the  tunnel  and  the  refrigeration  plant  operators  and  adjusted  to  match  as 
closely  as  possible. 

Test  section  static  (ambient)  temperature  was  obtained  by  subtracting 
from  the  measured  total  temperature  the  stagnation  temperature  Increment 


as  follows: 

*A  *  *1  *  26  "Cp  J 

(3-1) 
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true  airspeed,  ft/ sec 

2 

acceleration  of  gravity,  32.174  ft/sec 
mean  specific  heat  of  air,  0.24  Btu/(lb  -  °R) 
meachanical  equivalent  of  heat  in  engineering  units, 
778.2  ft-1 bs/Btu 

Liquid  water  content  (LUC)  capabilities  of  the  tunnel  are  limited  to 
some  extent  by  an  interdependence  on  airspeed  and  median  droplet  size 
because  higher  LWC  is  achievable  only  at  the  lower  airspeeds'.  Median 
droplet  size,  in  an  approximately  Langmuir  Vb"  size  distribution,  is 
somewhat  dependent  on  LWC  and  airspeed  in  that  larger  median  droplet 
sizes  cannot  be  obtained  at  a  combination  of  high  LWC  and  high  airspeed. 


Where :  V  = 
6  = 


J  = 


A  desired  combination  of  LWC  and  median  droplet  size  is  obtained  by 
setting  tunnel  controls  for  the  spray  system  air  and  water  pressures. 
To  determine  these  pressures  the  following  method  was  provided  by 
NASA: 


Given:  LWC,  0DM,  V  (e.g.,  V1  or  V2) 
Obtain: 

Ww  =  11.667  LWC(Vj) 

P  _P  . 

H20  air  38.7  \  77  / 


p  _P  .  j_  f 

H20  air  38.7  \  77  / 

,13[43.9-(4V2)-5]  C38.7(Ph  „-Pa1r)-S]  (  ^)'5 


(3-2) 


(3-3) 


bDM  -  0.006(VZ)  -  4 


-11.3  (3-4) 
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Where: 


LWC 

°DM 

V1 


■P. 


Ww 

3h2o 


air 

rair 


3 

-  liquid  water  content,  g/m 

*  median  droplet  size,  ym 

*  mph  true  airspeed 

=  knots  true  airspeed 

=  water  flow  rate,  Ib/hr 

2 

*  water  pressure,  Ib/in  gage 

2 

=  air  pressure,  lb/ in  gage 
=180  +  460  =  640°P 


The  tunnel  spray  system  water  and  air  pressures  were  calibrated  by 
NASA-Lewis  personnel. 
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3.2  Model  and  Instrumentation 

The  wind  tunnel  model  employed  for  icing  tests  was  adapted  from  a  full-scale 
nacelle  for  the  Garret-Ai Research  Corporation  TFE731-2  turbofan  engine. 
Modification  to  the  nacelle  included  incorporation  of  a  roll-formed  and 
welded  steel  duct  with  flanges  for  mounting  of  fore  and  aft  nacelle 
bodies  and  the  wrap  cowls.  A  strut  assembly  was  fabricated  from  square 
steel  tubes  which  was  then  welded  to  the  duct  and  to  mounting  plates  for 
installation  on  the  wind  tunnel  turntable.  The  nacelle  was  mounted  on  its 
side  in  the  tunnel,  as  shown  in  Figure  3-1 .  to  facilitate  angle-of -attack 
changes  by  rotation  of  the  turntable. 


Figure  3-1  Test  Nacelle  Installation 
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An  aluminum  plug  was  attached  to  the  nacelle  exit  for  certain  tests  to 
provide  a  variation  of  air  mars  flow  through  the  nacelle  and  consequently 
a  shift  of  the  stagnation  point  on  the  inlet  lip.  The  plug  was  adjustable 
in  fore  and  aft  position  through  use  of  threaded  attachment  rods.  The 
forward  cone  of  the  plug  was  electrically  heated  to  prevent  icing. 

Instrumentation  was  provided  on  the  nacelle  inlet  lip  to  measure  outside 
surface  temperatures  and  static  pressures  along  three  chordwise  positions* 
located  at  radials  of  6,  9,  and  12  o'clock.  Figure  3-2  illustrates 
typical  instrument  positions  on  the  lip.  These  positions  are  relative  to 
nacelle  orientation  as  installed  on  the  aircraft.  In  the  wind  tunnel, 
with  the  nacelle  mounted  on  its  side,  the  12  o'clock  position  faced  towards 
the  control  room  and  the  9  o'clock  position  was  on  the  top.  At  each 
radial  position  thermocouples  and  pressure  .aps  were  separated  circum¬ 
ferentially  by  approximately  two  (2)  inches.  Pressure  tubing  and  thermo¬ 
couple  wires  were  routed  inside  the  lip  chamber  in  order  to  present  a 
smooth  exterior  lip  surface. 

Hot  air  mass  flow  provided  to  the  inlet  lip  heat  distribution  system  was 
supplied  by  a  gas-heated  tunnel  supply  system.  Air  mass  flow  was  measured 
with  a  calibrated  orifice  located  outside  and  beneath  the  tunnel  test 
Section.  From  there  the  hot  air  was  ducted  up  through  the  test  section 
turntable.  Into  the  forebody  of  the  nacelle,  and  then  into  the  inlet  lip 
anti-icing  system.  Temperature  and  pressure  of  the  hot  air  entering  the 
anti-icing  system  were  measured  with  total  temperature  and  pressure  probes. 

Air  temperature  was  controlled  by  a  NASA  provided  Brown  Recorder/ 

Controller  which  regulated  a  natural  gas  fired  air  heater.  Air  pressure 
was  controlled  by  a  NASA  provided  pneumatic  servo  valve  which  adjusted 
a  remote  pressure  regulator. 


*  The  turbofan  engine  used  with  this  nacelle  produced  negllble  Inlet 
swirl. 
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All  pressures  and  temperatures  on  the  model,  tunnel  airspeed,  and  the 
heated  air  flow  provided  to  the  model  were  measured  and  recorded  with  the 
following  contractor  supplied  equipment: 

a.  Pressures  were  measured  using  a  scanning  valve  and  pressure 
transducers. 

b.  Temperatures  were  measured  using  iron-constantan  thermo¬ 
couples,  a  150°F  temperature  controlled  reference,  and 
amplifier. 

c.  Electrical  outputs  from  pressure  transducers  and  thermo¬ 
couples  were  measured  by  a  Hewlett  Packard  3455A  digital 
micro- voltmeter  with  an  accuracy  of  ±(0.007%  of  reading  +4 
digits)  in. the  0.1  volt  range. 

d.  Test  data  from  the  voltmeter  were  further  processed  with  a 
Hewlett-Packard  9825A  desk-top  digital  computer  to  apply 
instrument  calibrations,  convert  voltmeter  Outputs  to 
engineering  units,  calculate  tunnel  airspeed  and  hot  air 
mass  flow,  annotate  output,  record  data  on  magnetic 
cassette  tape,  and  print  outputs  on  paper  tape. 

Calibration  of  contractor  supplied  equipment  was  accomplished  prior  to 
the  Icing  tunnel  tests.  This  equipment  Included  all  Instrumentation 
employed  to  measure  model  temperatures,  pressures,  tunnel  air  speed, 
pressure  altitude,  and  hot  air  bleed  mass  flow  rate.  Tunnel  air  speed  was 
obtained  with  an  existing  pitot-static-probe  installed  In  the  test  section. 
Total  a^d  static  pressures  were  picked  off  the  tunnel  plumbing  system  for 
measurement  by  contractor  instrumentation. 

3.3  Scope  of  Tests 

The  scope  of  tests  conducted,  listed  In  Table  3-1,  were  performed  over  a 
period  of  six  weeks.  The  conditions  shown  in  Table  3-1  were  set  as  test 
objectives  and  were  not  always  what  was  achieved  during  this  test.  Actual 
test  conditions  are  presented  in  Appendix  C.  Installation  and  check-out 
of  the  test  model  and  instrumentation  required  four  days. 
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TABLE  3-1.  ICING  WIND  TUNNEL  TEST  SCHEDULE 
Dry  Air  Runs  -  Configuration  A* 

NOTE:  Conditions  shown  are  test  objectives.  Corrected  test  conditions 


are 

shown  in 

Appendix 

C. 

Run 

we 

°DM 

LWC 

o 

V 

Spray 

Time 

No. 

iiri 

{tb/niin) 

_USL 

MO. 

m 

(KIAS) 

(min.) 

1A 

3S0 

12 

-- 

— 

30 

225 

— 

2A 

350 

15 

— 

— 

30 

225 

-- 

3A  . 

400 

12 

— 

<*»  m 

30 

225 

m  m 

4A 

400 

15 

— 

— 

30 

225 

— 

5A 

450 

12 

— 

— 

30 

225 

— 

6A 

450 

15 

— 

— 

30 

225 

7A 

350 

12 

<m  m 

—  • 

14 

225 

mm 

8A 

350 

15 

mm 

mm 

14 

225 

mm 

9A 

400 

12 

14 

225 

— 

10A 

400 

15 

— 

•  •• 

14 

225 

mm 

I1A 

450 

12 

--  ' 

— 

14 

225 

mm 

12A 

450 

15 

— 

— 

14 

225 

mm 

13A 

350 

12 

m  m 

— 

-4 

150 

mm 

14A 

350 

15 

— 

-4 

150 

mm 

15A 

400 

12 

—  ' 

— 

-4 

150 

mm 

16A 

400 

15 

mm 

mm 

-4 

150 

'  mm 

17A 

450 

12 

— 

mm 

-4 

150 

mm 

18A 

450 

15 

mm 

— 

-4 

150 

mm 

19A 

350 

12 

•  «i 

w  V 

-15 

200 

mm 

20A 

350 

15 

mm 

•  «» 

-15 

200 

"im 

*  Single  skin  Inlet  lip  configuration. 


3-8 


Run 

lB 

WB 

LWC 

■5 

tT 

V 

Spray 

Time 

No. 

m 

(lb/min) 

um 

(g/m3) 

(°F) 

( KIAS) 

(min. ) 

21A 

400 

12 

-15 

200 

— 

22A 

400 

15 

-- 

-- 

-15 

200 

-- 

23A 

450 

12 

— 

— 

-15 

200 

-- 

24A 

450 

15 

.  — 

— 

-15 

200 

-- 

25A-28A  (Not 

used) 

29A**  400 

12 

— 

-- 

14 

225 

— 

30A**  400 

12 

— 

'  — 

14 

225 

— 

31At 

350 

12 

— 

— 

14 

225 

—  a  »  -2 

32At 

350 

12 

__ 

__ 

14 

225 

—  a  *  -4 

**  Variable  stagnation  point  using  nacelle  exit  plug, 
t  Variable  nacelle  angle  of  attack  conditions. 


ai  ion  n 

NOTE;  Conditions  shown  are  +*<■*■  L. 

are  shown  in  Appendix  C.  °  Jectives*  Corrected  test  condition: 


Run 

No^ 

*B 

.  m 

WB 

ilb/min) 

°DM 

LWC 

Isa si 

60A 

350 

12 

15 

0.78 

61A 

350 

15 

15 

0.78 

62A 

400 

12 

15 

0.78 

63A 

400 

15 

15 

0.78 

64A 

450 

12 

15 

0.78 

65A 

450 

15 

15 

0.78 

66A 

350 

12 

15 

0.6 

67A 

350 

15 

15 

0.6 

68A 

400 

12 

15 

0.6 

69A 

400 

15 

15 

0.6 

70A 

450 

12 

15 

0.6 

71A  450 

7?A  350 

73A  350 

74A  400 

75A  400 

76A  450 

77A  450 

.  78A  350 

79A  350 

80A  400 

81A  40Q 


m 

V 

(KXAS> 

Time 

iSlnJ. 

30 

225 

4.63 

30 

225 

4.63 

30 

225 

4.63 

30 

225 

4.63 

30 

225 

4.63 

30 

225 

4.63 

15 

12 

15 

12 

IS 

IT 

15 

12 

15 

12 

15 


15 


0.6 


14 

14 

14 

14 

14 

14 


225 

225 

225 

225 

225 

225 


4.63 

4.63 

4.63 

4.63 

4.63 

4.63 


20 

1.72 

-4 

150 

2.0 

20 

1.72 

-4 

150 

2.0 

20 

1.72 

-4 

150 

2.0 

20 

1.72 

-4 

150 

2.0 

20 

1.72 

-4 

150 

2.0 

20 

1.72 

-4 

150 

2.0 

20 

1.28 

-15 

200 

2.0 

20 

1.28 

-15 

200 

2.0 

20 

1.28 

-15 

200 

2.0 

20 

1.28 

-15 

200 

2.0 
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Wet  Air  Runs  -  Configuration  A  (Cont): 


Run 

No.  (°F) 

WB 

( 1 b/min) 

ddm 

ym 

LWC 

(g/m3) 

*T 

(°F) 

V 

(KIAS) 

Spr4y 

Time 
(min. ) 

82A  450 

12 

20 

1.28 

-15 

200 

2.0 

83A  450 

15 

20 

1.28 

-15 

200 

2.0 

84A-94A  (not  used) 

9 5 A*  400 

12 

15 

0.6 

14 

225 

4.63 

96A*  400 

12 

15 

0.6 

14 

225 

4.63 

97A**  350 

12 

15 

0.6 

14 

225 

4.63  a  =  -2° 

98A** , 350 

12 

15 

0.6 

14 

225 

4.63  a  *  -4 

*  Variable  stagnation 

point. 

**  Variable 
were  not 

nacelle  angle  of  attack, 
run  due  to  lack  of  time  in 

Note: 

tunnel 

These  two 
schedule. 

conditions 

NOTE:  Conditions ' shown  are  test  objectives.  Corrected  test  conditions 
are  shown  in  Appendix  C. 


Run 

No. 

h 

(°F) 

WB 

.  (Ib/min) 

ddm 

ym 

LWC 
( g/m3) 

*1 

(°F> 

V 

(K-IAS) 

358 

350 

12 

— 

— 

30 

225 

36B 

350 

15 

-- 

-- 

30 

225 

37B 

400 

12 

-- 

— 

.30 

225 

388 

400 

15 

-- 

30 

225 

39B 

450 

12 

-- 

— 

30 

225 

4  OB 

450 

15 

~ 

~ 

30 

225 

41B 

350 

12 

-- 

— 

14 

225 

42B 

.  350 

15 

-- 

14 

225 

43B 

400 

i 

..  w 

-- 

-- 

14 

225 

44B 

400 

15 

— 

—  . 

14 

225 

45B 

450 

12 

— 

-- 

14 

225 

46B 

450 

15 

— 

— 

14 

225 

473-52B  (not  used) 

53B 

350 

12 

—  1 

w 

-15 

200 

54B 

350 

15 

•mm 

— 

-15 

200 

55B 

400 

12 

— 

— 

-15 

200 

56B 

400 

15 

— 

-15 

200 

57B 

450 

12 

— 

mm 

-15 

200 

58B 

450 

15 

— 

— 

-15 

200 

*  Double  s kin  inlet  lip  configuration.  • 


Wet  Air  Runs  -  Configuration  B 

NOTE:  Conditions  shown  are  test  objectives.  Corrected  test  conditions 
are  shown  in  Appendix  C. 


Run 

No. 

h 

(°F) 

WB 

(lb/min) 

ddm 

ym 

LWC 

>  (g/m3) 

*1 

(°F) 

V 

(KIAS) 

Spray 
Time 
(min. ) 

100B 

350 

12 

15 

0.78 

30 

225 

4.63 

101B 

350 

15 

15 

0.78 

30 

225 

4.63 

102B 

400 

12 

15 

0.78 

30 

225 

4.63 

103B 

400 

15 

15 

0.78 

30 

225 

4.63 

104B 

450 

12 

15 

0.78 

30 

225 

4.63 

105B 

450 

15 

15 

0.78 

30 

225 

4.63 

106B 

350 

12 

15 

0.6 

14 

225 

4.63 

107B 

350 

15 

15 

0,6 

14 

225 

4.63 

108B 

400 

12 

15 

0.6 

14 

225 

4.63 

1093 

400 

15 

15 

0.6 

14 

225 

4.63 

HOB 

450 

12 

15 

0.6 

14 

225 

4.63 

111B 

450 

15 

15 

0.6 

14 

225 

4.63 

112B 

350 

12 

20 

1.28 

-15 

200 

2.0 

113B 

350 

15 

20 

1.20 

-15 

200 

2.0 

114B 

400 

12 

20 

l.:;8 

-15 

200 

2.0 

115B 

400 

15 

20 

1.28 

-15 

200 

2.0 

1166 

450 

12 

20 

1.28 

-15 

200 

2.0 

117B 

450 

15 

20 

1.28 

-15 

200 

2.0 

3-  13 


3.4  Icing  Wind  Tunnel  Test  Procedures 


Dry  air  testing  Was  conducted  prior  to  icing  tests  to  determine  the  inlet 
lip  skin  temperature  and  surface  static  pressure  profiles.  Various  combi¬ 
nations  of  tunnel  airspeed,  tunnel  temperature,  hot  air  bleed  mass -flow  rate 
and  temperature,  and  nacelle  angle-of-attack,  as  indicated  in  Table  3-1, 
were  run  with,  dry  air  conditions.  Following  the  dry  runs,  wet-air  tests 
were  conducted  for  the  same  conditions  except  that  cloud  characteristics 
were  simulated  using  the  tunnel  water  spray  system. 

The  procedure  employed  during  dry  air  runs  is  described  as  follows: 

a.  Start  tunnel  and  set  at  idle  speed. 

b.  Start  tunnel  cool-down. 

c.  When  tunnel  has  reached  the  approximate  temperature  desired, 
increase  airspeed  to  test  condition  and  re-adjust  tunnel  temp- 
ature. 

d.  While  (c)  is  in  progress,  fire  air-bleed  gas  heater. 

e.  Adjust  air  heater  temperature  and  bleed-air  pressure  to  obtain 
desired  mass  flow  rate  and  temperature  at  entrance  to  nacelle 
inlet  lip  anti-icing  system. 

f.  When  all  conditions  are  stabilized,  conduct  instrumentation  data 
scans,  convert  to  engineering  units,  and  record  data. 

g.  Hold  tunnel  conditions  constant  and  adjust  bleed  air  temperature 
or  mass  flow  rate  to  next  test  point  and  again  stabilize  and 
record  data. 

h.  When  all  test  conditions  at  a  given  tunnel  airspeed  or  temperature 
are  complete,  adjust  to  next  test  point  and  repeat  procedure  above. 

Following  completion  of  dry  runs  for  a  particular  test  model  configuration 
the  procedure  is  repeated  for  wet  runs,  except  that  the  tunnel  water/air 
spray  system  is  used  to  simulate  icing  cloud  conditions.  Dry  air  test 
procedures  are  repeated  through  item  (g)  with  the  following  additional  steps 
used  during  icing  test  runs:  .  * 


h.  When  all  conditions  are  stabilized,  including  p re -adjustment  of 
spray  system  water  and  air  pressure,  switch  spray  system  on, 
adjust  water  flow  rate  rotameters  for  uniform  spray  across  tunnel 
cross  section,  record  data  after  a  predetermined  interval,  and  at 
completion  of  spray  time  switch  spray  system  and  hot  air  bleed 
flow  off. 

i.  As  soon  as  spray  system  and  bleed  air  are  off,  conduct  a  rapid 
shut-down  of  the  tunnel  fan  motor.  This  procedure  freezes  the 
"iced"  condition  of  the  model  because  tunnel  temperature  has  a 
tendancy  to  cool-down  further  by  10-15°F  and  to  hold  this 
condition  due  to  the  cold-soaked  state  of  the  all  metal  tunnel. 

j.  Conduct  an  inspection  of  the  ice  formation  on  the  model,  take 
measurements  and  samples  if  desired,  and  remove  ice  from  the  model 
using  a  steam  hose. 

k.  Following  the  above,  subsequent  runs  are  conducted  repeating  the 
procedure  outlined. 

Spray  duration  during  a  test  run  was  selected  corresponding  to  the  approxi¬ 
mate  time  it  would  take  to  fly  through  the  cloud  size  being  simulated  at 
the  tunnel  airspeed.  During  simulation  of  maximum  continuous  cloud 
characteristics,  corresponding  to  low  levels  of  LWC,  tunnel  airspeed  was 
run  at  225  KIAS  (approximately  230-235  KTAS,  depending  on  tunnel 
temperature). 

For  intermittent  maximum  cloud  simulation  conditions,  however,  airspeed  was 
limited  to  either  200  KIAS  (approximately  196  KTAS)  or  150  KIAS  (approxi¬ 
mately  142  KTAS)  because  the  tunnel  spray  system  is  limited  in  maximum 
water  flow  rate.  Run  time  based  on  the  time  of  transit  through  intermittent 
cloud  size  at  these  speeds,  however,  is  insufficient  to  allow  data  acquisition 
so  run  time  was  increased  to  two  (2)  minutes  as  was  seen  In  Table  3-1. 
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Although  time  in  a  cloud  may  have  no  significance  in  developing  an 
analytical  method  there  is  other  rationale  for  limiting  the  exposure  time. 
For  the  nacelle  and  anti-icing  system  employed  there  are  design  limits  to 
which  icing  can  be  prevented  at  any  bleed  air  temperature  and/or  mass 
flow  rate  when  testing  in  conditions  of  either  maximum  continuous  or  maximum 
intermittent  icing.  It  was  to  avoid  exceeding  these  limits  that  the 
approach  was  taken  to  use  FAA  defined  cloud  sizes  and  tunnel  velocity  to 
determine  exposure  time.  In  addition,  longer  exposure  time  to  the  icing 
cloud  results  in  an  undesirable  amount  of  ice  forming  on  the  unheated 
nacelle  support  structure,  especially  under  glaze  icing  conditions. 


4.0  CORRELATION  OF  RESULTS 


Dry  air  wind  tunnel  test  data  were  used  with  the  computer  program  CHANEFF 
to  determine  effective  channel  efficiencies  for  each  set  of  data  corres¬ 
ponding  to  a  different  nacelle  lip  pressure  coefficient  distribution  and 
internal  heat  distribution  system.  These  are  related  to  changes  in  nacelle 
angle  of  attack,  nacelle  air  flow,  circumferential  position  around  the 
lip  and  double  versus  single  skin  internal  distribution  systems.  Although 
it  would  appear  that  a  change  of  internal  bleed  air  flow  or  bleed  air 
temperature  would  result  in  different  channel  efficiencies,  this  is  not 
the  case.  For  a  given  type  of  internal  distribution  system,  the  surface 
temperature,  which  is  used  in  determining  the  channel  efficiency,  is 
directly  related  to  the  quantity  of  heat  provided  whether  by  virtue  of 
higher  bleed  mass  flow  or  bleed  air  temperature.  Thus  the  calculated 
efficiency,  for  the  range  of  bleed  flows  and  temperatures  tested,  is 
independent  of  these  two  parameters.  Hence  in  order  to  use  a  larger  data 
base  for  improved  accuracy  of  the  channel  efficiency,  the  individual  test 
efficiencies  for  a  given  physical  arrangement  are  averaged.  There  will 
be  efficiencies  related  to  the  6,  9  and  12  o'clock  measurement  positions 
for  each  variation  in  nacelle  angle  of  attack,  mass  flow  and  internal 
distribution  system.  In  total  18  different  efficiency  distributions  were 
computed,  corresponding  to  three  angles  of  attack,  two  mass  flows  and  two 
distribution  systems  for  each  of  three  clockwise  measuring  stations.  Table 
4-1  presents  a  matrix  of  which  dry  test  runs  were  used  to  calculate  the 
efficiency  and  to  which  wet  test  runs  they  apply. 


Table  4-1  Efficiency  -  Run  matrix 

■ 

Runs  used 

Runs  EFF  applies  for  T$ 

Calculations 

for  EFF 

Dry 

Wet 

10A-24A 

1A-24A 

60A-83A 

29A 

29A 

95A 

30A 

30A 

96A 

31A 

31A 

- 

32A 

32A 

- 

35B-4SB 

35B-46B 

• 

CO 

8 

iH 

53B-58B 

53B-58B 
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A  typical  averaged  surface  pressure  coefficient  is  presented  in  Figure  4-1. 

It  is  noted  that  the  peak  pressure  coefficient  from  this  test  data  is 
greater  than  1.0.  Since  this  is  physically  impossible  the  Pc  curve  was 
faired  to  1.0  at  the  stagnation  point  as  shown.  This  appeared  to  be  a 
common  occurance  with  the  measured  pressure  data.  This  first  came  to 
attention  during  the  running  of  the  test  and  an  effort  was  made  to 
determine  the  cause  of  it.  No  reason  for  its  occurance  would  be  determined. 
Corresponding  measured  surface  temperature  distributions  are  shown  in 
Figure  4-2  and  the  resulting  averaged  channel  efficiencies  from  program 
CHANEFF  are  presented  in  Figure  4-3.  Utilizing  the  pressure  data  from 
Figure  4-1  and  efficiencies  from  Figure  4-3  surface  temperatures  were 
computed  with  program  ICEOFF.  A  comparison  of  calculated  and  measured 
surface  temperature  for  these  dry  air  conditions  is  presented  in  Figure 
4-4.  It  is  seen  that  the  agreement  is  excellant. 

A  comparison  of  surface  pressure  coefficient  at  the  three  nacelle  inlet 
lip  measuring  stations  is  shown  in  Figure  4-5  for  both  wet  and  dry  air.  As 
would  be  expected  the  presence  of  water  has  little  if  any  effect  on  the 
pressure  coefficient.  The  corresponding  channel  efficiencies  for  dry  air 
for  the  three  inlet  lip  measuring  stations  is  given  in  Figure  4-6  and  the 
comparison  of  calculated  and  measured  surface  temperature  for  the  three 
inlet  lip  measuring  stations  for  both  single  and  double  skin  bleed  air  distri¬ 
bution  systems  is  presented  in  Figure  4-7.  It  Is  noted  in  Figure  :  -6 
that  at  the  6  o'clock  position,  the  calculated  channel  efficiency  t.xceeds 
1.0.  Although  it  is  physically  impossible  to  exceed  190*  efficiency,  the 
mechanics  of  the  calculation  process  allow  this  to  happen  on  rare  occasions. 
Detailed  checks  revealed  no  mathematical  or  computational  error  and  the  use 
of  this  efficiency  in  program  ICEOFF  gave  an  exart  reproduction  of  the 
measured  surface  temperature.  It  is  concluded  that  there  must  be  some  of 
the  basic  assumptions  underlying  the  method  of  channel  efficiency  that 
allows  this  to  happen.  Examination  of  the  equation  used  to  compute  EFF 
reveals  how  this  can  occur.  From  equation  (2-48),  (2-57),  (2-58)  and 
(2-60): 
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It  can  be  shown  that,  compared  to  Tg  and  T^,  T ^  and  Tg  are  relatively 
unimportant  for  the  test  conditions  of  this  report.  Hence  it  is  easily 
seen  that  for  conditions  such  that  Tg  and  Tg  are  quite  close  in  magnitude 
and  Tg  and  are  quite  different  in  magnitude,  the  result  for  EFF  could 
exceed  1.0.  This  is  undoubtedly  due  to  some  of  the  simpl ifications  that 
resulted  from  early  assumptions  concerning  the  flow  and  heat  transfer 
characteristics.  Understanding  how  it  can  occur  makes  it  no  more  desirable 
but  in  the  final  analysis  it  will  reproduce  accurate  surface  temperatures 
and  thus  will  be  tolerated.  Examination  of  Figures  4-4  and  4-7,  reveals 
that  the  agreement  for  dry  air  is  excellent  but  for  wet  air  it  is  less 
good.  It  appears  that  the  method  tends  to  over  predict  the  temperature 
near  the  stagnation  point  and  to  under  predict  It  near  the  aft  limit  of 
the  heated  area.  Although  the  amount  of  difference  between  calculated 
and  measured  Tg  varies  with  various  conditions  the  trend  noted  above  is 
fairly  consistant.  This  could  be  caused  by  a  somewhat  inaccurate  modeling 
of  either  the  evaporation  process  or  the  heat  transfer  process  that 
resulted  from  simplifying  assumptions. 


The  balance  of  the  figures  in  this  section  present  the  effects  of 
independent  parameters  on  the  pressure  coefficient,  the  resulting  effic¬ 
iency  and  a  comparison  of  calculated  and  measured  surface  temperature. 
Figures  4-8,  4-9  and  4-10  present  this  information  as  influenced  by 
nacelle  angle  of  attack.  Only  dry  air  temperatures  are  presented  because 
time  did  not  allow  the  testing  of  the  wet  air  counterparts.  ' 

Figures  4-11,  4-12  and  4-13  present  pressure  coefficient,  efficiency  and 
temperature  comparison  for  changes  In  inlet  mass  flow.  Here  again  one  of 
the  efficiency  curves  (Run  29A)  exhibits  a  value  greater  than  1.0.  In 
this  case  2.44.  Although  appearing  to  be  unreasonable  at  first,  the  same 
arguments  as  set  forth  before  apply  and  again  the  calculated  dry  air 
surface  temperature  agrees  well  with  measured  values. 
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Figures  4-14  and  4-15  present  channel  efficiency  and  comparison  of  temp¬ 
eratures  for  single  and  double  skin  anti -icing  systems. 

The  last  four  figures  (4-16,  4-17,  4-18  and  4-19)  show  only  comparative 
temperatures  because  they  were  aerodynamically  and  geometrically  similar 
to  other  configurations  and  used  pressure  coefficients  and  efficiencies 
defined  for  them.  These  similarities  are  distinguished  in  fable  4-1. 

The  effects  of  bleed  air  temperature  are  illustrated  in  Figure  4-16.  Here 
the  comparison  include:  wet  and  dry  air  as  well  as  single  and  double  skin 
systems. 

Temperature  comparison  for  changes  in  bleed  air  mass  flow  are  shown  in 
Figure  4-17.  Included  are  data  for  wet  and  dry  air  as  well  as  single  and 
double  skin  systems. 

Figure  4-18  gives  comparative  temperature  for  a  variety  of  liquid  water 
contents.  It  should  be  noted  here  that  it  was  not  possible  to  hold  all 
other  test  parameters  constant  during  the  test  while  varying  liquid  water 
content  alone.  Hence  a  direct  comparison  for  the  effects  of  liquid  water 
content  alone  cannot  be  made.  Again  both  single  and  double  skin  systems 
are  shown. 

The  comparison  of  temperature  while  changing  water  droplet  diameter  Is 
given  In  Figure  4-19  for  both  single  and  double  skin  systems.  As  for 
liquid  water  content.  It  was  not  possible  to  change  droplet  diameter  alone 
while  holding  all  other  test  parameters  constant.  Thus  a  direct  comparison 
to  determine  effects  of  droplet  diameter  alone  Is  not  possible. 

For  all  of  these  comparisons,  the  results  are  nearly  the  same:  very  good 
agreement  for  dry  air  and  moderately  good  agreement  for  wet  air.  It  Is 
noted  that  even  minute  as  well  as  large  fluctuations  In  surface  temperature 
are  reproduced  for  dry  air  calculations.  Although  not  as  accurate  as  fop 
dry  air,  the  wet  air  temperature  distributions  generally  give  the  correct 
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trend  and  do  display  the  gross  variations  that  exist  in  the  measured  data. 
This  leads  to  the  belief  that  the  method  is  very  useful  and  needs  only 
some  refinements  as  applied  to  wet  air  calculations. 

A  correlation  of  all  calculated  and  measured  data  was  performed.  For  the 
966  calculated  dry  air  temperatures  the  mean  deviation  from  measured 
temperature  is  6°F  and  for  the  924  wet  air  temperatures  the  mean  deviation 
is  19°F. 

During  the  data  analysis,  it  was  noted  that  the  calculated  surface 
temperatures  for  wet  air  consistantly  gave  better  correlation  with  the 
measured  results  for  the  double  skin  system  than  for  the  single  skin 
system.  There  are  at  least  two  reasons  why  this  is  so.  First  the  double 
skin  will  give  a  more  uniform  temperature  distribution  (i.e.  nearly 
constant)  and  second,  the  deviations  from  this  nearly  constant  value  will 
be  smaller.  Both  of  these  will  work  in  favor  of  better  accuracy  from  the 
calculation  method.  These  are  well  illustrated  by  the  tabulated  data, 
(Tables  C-4,  C-5  and  C-6)  and  Figures  4-17  and  4-18. 
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Figure  4-1  (Cent.) 


4-7 


SURFACE  TEMPERRTURE  -  L  -  DEG  F  SURFACE  TEMPERATURE  *  \  -  DEG  F 


DISTRNCE  ALONG  SURFACE  -  INCHES 
Figure  4-2  Typical  Measured  Surface  Temperature  Distribution 
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Figure  4-3  Average  Channel  Efficiency 
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Figure  4-4  Comparison  of  Calculated  and  Measured  Surface  Temperature 
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Figure  4-6  Effects  of  Nacelle  Lip  Position  on  Efficiency  Distribution 
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Figure  4-7  Comparison  of  Calculated  and  Measured 
Surface  Temperature  for  Three  Nacelle 
Llo  Positions. 
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Figure  4-7  (Cont.) 
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Figure  4-7  (Cont.) 
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Figure  4-8  Effects  of  Nacelle  Angle  of  Attack 
on  Surface  Pressure  Coefficient 
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Figure  4-9  Effects  of  Nacelle  Angle  of  Attack 
on  Channel  Efficiency 
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Figure  4-10  Comparison  of  Calculated  and  Measured 
Surface  Temperature  for  3  Nacelle 
Angles  of  attack 
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Figure  4-11  Effects  of  Nacelle  Inlet  Mass  Flow 
on  Surface  Pressure  Coefficient 
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Figure  4-13  (Cont.) 
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Figure  4-17  Comparison  of  Calculated  and  Measured 
Surface  Temperature  for  Two  Bleed 
Air  Flow  Rates. 
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Figure  4-17  (Cont.) 
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Figure  4-18  Comparison  of  Calculated  and  Measured 

Surface  Temperature  for  Different  Liquid 
Water  Content. 
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Figure  4-18  (Cont.) 
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Figure  4-19  Comparison  of  Calculated  and  Measured 
Surface  Temperature  for  Different  Water 
Droplet  Diameters.  -  Single  skin. 
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Figure  4-19  -  (Cont.)  -  Double  Skin 


5.0  APPLICATION  OF  THE  PROCEDURES 

The  analysis  method' consists  of  two  separate  programs  written  in  Fortran 
IV.  As  presently  programmed,  these  can  be  run  on  PDP  11/70,  IBM  370  or 
other  compatible  equipment. 

The  first  program,  CHANEFF,  is  used  to  determine  the  heat  transfer 
efficiency  of  the  inlet  lip  hot  air  distribution  system.  The  intent  is 
to  determine  the  efficiency  for  the  particular  configuration  to  be 
employed  in  the  nacelle  design  by  running  dry  air  tests  in  a  wind  tunnel. 

The  test'  model  can  be  two-dimensional  with  a  configuration  representing  a 
partial  section  of  the  nacelle  inlet  design.  .Data  to  be  collected  during 
dry  air  tests  and  geometrical  information  that  is  required  for  running 
program  Ch.  ,_FF  are:  . 

Nacelle  Geometry 

Leadii  <  edge  radius 

Streamwise  surface  length  of  the  heated  surface 
High-light  diameter 
Flight  Condition  Data 

Flight/tunnel  pressure  altitude 
Ambient  temperature 
Flight  velocity/airspeed 
Bleed  air  mass  flow  rate 
Bleed  air  temperature 
Inlet  Lip  Data 

Surface  skin  temperature  profile  streamwise  around  inlet  lip 
Surface  static  pressure  profile  streamwise  around  inlet  lip 

Surface  static  pressures  will  be  used  to  calculate  the  pressure  coefficient 
characteristics  of  the  inlet  lip.  With  the  above  listed  data  the  channel 
efficiency  can  be  calculated  for  incremental  intervals  streamwise  around 
the  inlet  lip.  For  the  same  nacelle  geometry  and  engine  flow  characteristics, 
the  pressure-coefficient  profile  will  remain  the  same  so  long  as  the  nacelle 
angle-of-attack  is  held  constant. 


The  second  program,  ICEOFF,  can  be  employed  to  evaluate  the  anti-icing 
performance  of  the  nacelle  inlet.  This  program  can  be  used  in  any  one  or 
all  of  the  following  ways: 

a.  As  a  tool  during  the  design  phase  to  determine  the  bleed  air 
requirements  for  anti-icing. 

b.  With  incorporation  of  an  engine  performance  routine  into  the 
present  program  this  can  be  used  to  size  the  heating  system. 

c.  As  a  means  to  demonstrate  engine/nacelle  anti-icing  capability 
as  part  of  the  FAR  anti-icing  performance. 

The  data  required  for  raining  program  ICEOFF  are: 

Nacelle  and  Aircraft  Geometry 
Nacelle  chord  length 
Nacelle  thickness-to-chord  ratio 
Nacelle  inlet  leading  edge  radius 
Nacelle  inlet  high-light  diameter 
Extent  of  inlet  lip  heated  surface 
Wing  area. 

Flight  Conditions 

Flight/ tunnel  pressure  altitude 
Aircraft  gross  weight 
Flight  velocity 

Constants  of  the  airplane  drag  polar 
Inlet  Lip  Data 

Surface  pressure  coefficient  profile  streamwise  around  the  Inlet  lip 
Channel  efficiency  profile  streamwise  around  the  inlet  lip  (From 
Program  CHANEFF) 

Icing  Conditions 

Ambient  temperature 
Droplet  median  diameter 
Liquid  water  content 
Cloud  horizontal  extent 


Program  ICEOFF  is  written  to  operate  in  two  different  modes  as  follows: 
Mode  I  -  This  mode  uses  icing  tunnel  data  and  predicts  evaporation 
and  runback.  In  order  to  check  run-data  as  it  proceeds,  this  mode 
of  operation  prints  out  intermediate  calculation  of  the  point-by¬ 
point  analysis. 

Mode  II  -  This  mode  is  used  for  icing  prediction  under  given 
flight  and  icing  conditions. 

Detailed  instructions  for  use  of  these  programs  along  with  sample  inputs 
and  outputs  are  provided  in  Appendices  A  and  8. 
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6.0  CONCLUSIONS 


Two  computer  programs,  based  on  a  simpl if ied  heat  transfer  theory,  have  . 
been  developed.  The  channel  efficiency  program,  CHANEFF,  provides  a  means 
for  evaluating  the  external  and  internal  heat  transfer  characteristics  of 
an  engine  inlet  anti-icing  system.  The  anti-icing  program,  ICEOFF,  provide 
a  means  of  evaluating  the  icing  performance  characteristics  of  an  engine 
inlet  anti -icing  system.  Both  programs  are  developed  for  a  continuous 
flow  hot  gas  anti-icing  system. 

The  program  ICEOFF  calculates  the  inlet  lip  skin  temperature  profile  for 
dry  air  with  very  good  agreement.  The  calculations  for  icing  conditions, 
though  useful,  are  not  so  good.  The  mean  deviation  of  the  calculated 
temperature  is  6°F  for  dry  air  and  19°F  for  wet. 

It  is  concluded  that,  although  better  accuracy  could  bo  desirable  for  wet 
air  calculations,  the  method  is  simple  to  use  and  gives  rapid  and  useable 
predictions  for  the  surface  temperature  distribution  on  an  engine  nacelle 
inlet  lip  that  is  anti-iced  with  hot  bleed  air  flow.  The  method  is  applic¬ 
able  to  a  wide  range  of  independent  variables  in  the  problem  and  treats 
double  skin  distribution  systems  equally  as  well  as  single  skin  piccolo 
tube  systems. 
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7.0  RECOMMENDATIONS 


The  .underlying  basis  of  the  present  approach  is  the  characterization 
of  the  internal  heat  transfer  that  is  used  in  the  analysis  published 
by  the  Federal  Aviation  Administration1,  by  a  single  parameter,,  referred 
to  as  Channel  Efficiency.  This  parameter  is  used  in  lieu  of  a  more 
definite  determination  of  an  internal  convective  heat  transfer  coefficient 
Channel  Efficiency  equates  to  the  ratio  of  external  heat  transfer  divided 
by  internal  heat  transfer.  Any  errors  caused  by  improper  determination 
of  the  internal  or  external  heat  transfer  will  then  appear  in  the  value 
of  Channel  Efficiency  when  using  the  program  for  evaluating  experimental 
data. 

Based  on  the  analysis  of  dry  air  data  using  program  ICEOFF,  it  is 
concluded  that  the  assumptions  made  in  formulating  the  internal  and 
external  flow  and  heat  transfer  model  are  very  reasonable.  This  is 
supported  by  the  6°F  mean  temperature  deviation  which  is  probably  well 
within  the  data  accuracy.  This  is  considered  excel lant  agreement  and 
no  change  to  this  portion  of  the  program  is  recommended. 

The  correlation  accuracy  for  wet  air  (19°r  mean  deviation)  is  not  as 
good  as  for  dry  air.  It  appears  that  there  is  possibility  for  improve¬ 
ment  here.  Brief  analyses  of  the  cause  of  major  temperature  differences 
between  calculateu  and  measured  near  the  stagnation  point  indicated  that 
the  rate  of  heat  transfer  in  this  region  is  consistantly  over  predicted 
resulting  in  higher  than  measured  temperature.  On  the  other  hand,  near 
the  aft  limit  of  the  heated  region,  the  calculated  temperatures  tend  to 
be  lower  than  measured.  A  more  in-depth  investigation  of  these  two 
regions  would  probably  result  in  Improved  calculation  procedures  and 
better  accuracy  on  temperature  correlation.  Potential  candidates  for 
improvement  would  be  the  modeling  for  evaporation  and  heat  transfer. 


1.  op.  cit. 
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The  method  presented  does  give  accurate  trends  for  a  variety  of  external 
ambient  conditions  (i.e.  droplet  diameter,  liquid  water  content,  etc.) 
and  no  further  improvements  are  recommended  in  this  area. 
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This  program  uses  experimental  data  for  flight  and  bleed  air  conditions  to 
compute  a  channel  heat  transfer  efficiency  on  a  point  by  point  analysis. 

INPUT  PREPARATION 

CARD  1  (40A2)  Title 

This  card  can  contain  up  to  80  columns  of 
alphanumeric  information.  It  will  be  used  as 
the  headin  for  the  output  data  sheet. 

CARD  2  (4F8)  Nacelle  Ceometry 

This  card  has  4  fields  of  8  columns  each  to 
contain  the  variables  RLE,  SH,  DIA  and  DS 
(F8.2  format). 

Column  1  RLE  -  Radius  of  curvature  for  the  leading  edge  of 

the  nacelle  inlet  lip  in  feet. 

Column  9  SH  Streamwise  length  of  the  heated  surface  in 

feet.  (Heated  area  per  foot  circumference). 

Column  17  DIA  -  Nacelle  inlet  highlight  diameter  in  feet. 

Coliann  25  DS  -  Length  of  increment  used  between  points  of 

analysis  in  feet. 

CARD  3  (16F5)  Pressure  Coefficient  Table. 

This  card  has  16  fields  of  5  columns  each  to 
contain  streamwise  distance  measurements, 
from  the  stagnation  point  to  the  point  under 
analysis.  In  Inches.  The  values  should  begin 
at  zero  (stagnation  point)  and  Increase  at 
any  desired  Interval  until  the  aftmost  heated 
region  Is  reached  (F5. 2  format). 


A.  2  Input  and  Operation  (Cont.) 


CARD  4  ( 16F5) 


CARD  5  (40A2) 


CARD  6  (5F8) 


Col umn  1 
Column  9 
Column  17 
Column  25 
Column  33 

CARD  7  (16F5) 


CARD  8  (16F5) 


Pressure  Coefficient  Table. 

This  card  has  16  fields  of  5  columns  each 
to  contain  pressure  coefficient  values  for  the 
corresponding  points  on  Card  3  (F5.2  format). 

Flight  Title 

Test  or  flight  nunber  and  any  pertinent 
information  up  to  80  columns  (40A2  format). 

Flight  Parameters. 

This  card  has  5  fields  of  8  columns  each  to 
contain  the  variables  ALT.  TAP,  VKTAS, 

WBM  and  TBF  (F8.2  format). 

ALT  -  Altitude  (MSL)  of  aircraft  in  feet. 

VKTAS  -  True  airspeed  of  the  aircraft  in  knots. 

TAF  -  Ambient  temperature  in  BF. 

WBM  -  Bleed  air  mass  flow  per  engine  in  lb/min. 

TBF  -  Bleed  air  temperature  in  °F. 

Surface  Temperature  Table 
This  card  has  16  fields  of  5  columns  each  to 
contain  streamwise  distance  measurements,  from 
the  stagnation  point  to  the  point  under  analysis 
in  inches.  The  values  should  begin  at  zero 
(stagnation  point)  and  increase  at  any  desired 
interval  until  the  aftmost  heated  region  is 
reached  (F5.2  format). 

Surface  Temperature  Table. 

This  card  has  16  fields  of  5  columns  each  to 
contain  surface  temperature  values  corresponding 
to  the  points  on  Card  7  (F5.1  format). 


Up  to  100  cases  of  cards  5,  6,  7  and  8  can  be  handled  at  one  submission. 

Card  8  followed  by  a  blank  card  will  cause  the  routine  to  calculate  point- 
by-point  average  efficiencies  and  then  loop  back  to  read  new  nacelle 
geometry  and  pressure  distribution  data  (cards  2,  3  and  4).  As  many  as 
20  such  configurations  can  be  handled  at  one  submission.  Card  8  followed  by 
2  blank  cards  will  terminate  the  program. 


non 


t.3  Program  Listing 


t-BuDRA*  CHANF.FF . o/2>/7v . . . . 

t^GIr.f  IMc.1  I  IP  JMl-ICt  AI“AI  *SiS.... . 

L*  f  i  F  R  '*  I  N  AT  1  Oil  OF  F*LFtr  AIR  C  n  *  f«  h  fr.  L  FFKIC  IK  :<Ci 


RtAL,  J.Kl.Mj 

DlftKSTON  T I  ILK  (401  ,  F  L.  T  (  4  0  ) 

DI^f.MSION  P.b.DAl  (  1b,2) 

DIMtVSIOM  Tf>UAT(16,2> 

DI’JfN6lOK  F  K  F  1 C  (  1  b  ) 

DrRfcNSTON  M'CR(lo) 

OFF  Mf  UNITS?  .  UrtMfc.=  '  CkiANEF  F  .DAT  '  ,1  VFF  =  '01  i> '  ) 
Opt-  W  (  0N1T  =  3  -,NA*"e.=  '  CDA,'>'F  F  F  .OUT  •  ,L/IS?=  '  PRJ  Ml  •  ) 

”-r(i 'i  of  «kt 

TK  =  2 
Tw=3 
NP.Ss)  r 
NTS=1 b 
G=32 .  174 
J=77R. 

Pl=3.1416 

R=53.3S 

- - - TuF  UTS 


PEAOf JR.400]  T I  Hit 
DO  170  Ji.lsl  ,  2u 
X -i  =  0.0 
DO  b  U=l,7 
FFF]C(L)=0.0 
5  Cunt i  n  ii  t 

PEM>(1R,300)  FLt.iH.DI  A  ,D.S 
IFCftl.t  .It.  0.0)  CALL  EX JT 
SL-SH / 2 . 

NsSL/CS+1 


READ ( IP , 1 400 )  PbDAT 
DO  150  11=1.100 
READ(  IK, 400)  FI1 

READ! IP  ,  300)  A  L  T  ,  1  A  F’ ,  V  M  A  5  .  WBM ,  TBF 
IKCVKTAS  .It.  0.0)  (,fi  TO  160 
JFCVKfAS  .fit.  1.0)  GO  TO  20 
XMsVKTAS 
20  CONTINUE 
NCA5E=I1 


READ! IB, 1400)  TSOaT 
DO  25  L=  1  ,  1.6 

TSDAT(L,?)='iSUAHL,2)>45R.6B8 
25  CONTINUE 

TAsTAF+459.b«H 
TB=TBE+459 . feBb 
Wd=WBM/(60. *PI*DIA) 
teBlTF.(I»,4lO)  TITLE 
teKTTF(Ite,600)  RLE  ,  SKi  ,  Dl  A 
WHIT Elite, 400)  ELT 

teRIT Elite, 700)  ALT,TAF,VtTAS,WBM,T3F 
teKIT Elite, bOO) 


- - - r— CALCULATION  DF  BASIC  PAFtAF.ETFRS 

AsSOR! <1.4*G*R*1A) 

IKX*  .LE.  0.0)  Gu  TO  30 
VKT  AS  =  Xf*  *  A 
30  VQ*VKlAS*l.b87« 

CPs. 2365+7. b*TA/10.**6 
MUs7,475*TA**l ,5/(TA+21b. )/10.**7 
Kls(4, 72?+. 06944* TA)/10.**7 
PR=CP*MU/K1 

,  PA  =  2U6,21/t  XK  ALT/ (?77 10. -.098774*  ALT)  ) 
RHOsPA/TA/P 
REFTn=RHO*VO/*U 
P£LE*RFf TO*2.*PLF 
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A . 3  Program  Listing  (Cont.) 

dp=rnu/2 .  ♦vonn/o 
T2  =  vn*vn*HR**.  5/ ( 2.*cp*g*J) 

- -  - POI^T-RY-pntNT  ANALYSIS  OF  hF  AT  ThANSFF.R 

RELATIONSHIPS,  STARTING  A 1  STAGNATION 
PuTnT  AND  KNOJNG  a  aFTMOST  PCiNT  OF 
HEATED  REGION 

AT  AFTMOST  POINT  OF  HE»TFP  REGION 


- INITIALISE  CONDITIONS  A  N I)  S'i  API 

SsPSDATI 1 , 1 ) /l? 

N  N  =  1 

DO  1 5  0  Isl.u 

IF  (J  .GT.  1)  GO  TO  35 

- - —  ----  —  STAGNATION  POINT  CALCULA'llOf  S 


PC=RSDAT(  I  ,2) 

TS=TSD AT (1,2) 

THFTA=S/PLF*57. JO 
WRITFC  IW,90uKS,pr,  ThFTA 
PL=PA+DP 
CO  TO  40 
35  CONTINUE. 

- CALCULATIONS  AWAY  FRO,;  THE  ST AGNATION  PiJli.T 

PC  =  TRP(S.-PSL»AT,NPS) 

TSaTRPCS.TSDAT .NTS) 

THETA=S/RLFK57.30 
WRITE  ( IV>  ,yOO)S,PC,THFJTA 
V=VO*SORT( 1 .-PC) 

PL=PA+PC*DP 
RF.SsPEFT0*S*V/V0 
40  CONTINUE 

TF  (THETA  .GT.  25.0)  GO  TO  50 

—  — - - - — CYLINUFR  REGION 

Hs,57*Kl*pR»*.4*RELE**.5/HLfc*( I (THK1 A/ 90. ) **3 ) 
WRITE(IW,4?0) 

GO  TO  105 
50  CONTINUE 

IF( RES  . G'l  .  2.0L5)  GO  TO  60 

...... - -  —  —  laminar  FLOW  REGION 

Hs. 332*K1*PR**(  1 ./ 3. )*R£S**,5/S 
WRI1F(IW,430) 

GO  TC  105 
60  CONTINUE 

IF ( RFS  .GT.  1,?*J6)  GU  TO  100 
- - TRANSITION  FLOW  REGION 


S1 1  =S 

REST  1 1 *RES 
S1=S11-0S 
GO  TO  (70,95) ,  NN 
70  CONTINUE 


.............. CALCULATION  OF  PC,  V  AND  H  AT  THE  BEGINNING  OF 

THE  TRANSITION  REGION 

PC*TRP ( S 1 .PSDAT.NPS) 

VsVO*SQRT(l»-PC) 

REST1*RLFT0<IS1*V/V0 

IF(ABS(REST1-2,UF5)  .LF.  100.)  GO  Tl.  75 
SlIsSll  +  (Sl-Sll  )/(HF.STl-RfcSTii)*(2.0fi.5-RESTin 
SllaSl 

REST1 1*RFST1 
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A . 3  Program  Listing  (Cont.) 


51  =6  II 

C-0  10  7  0 
7  5  CONTINUE 

HLAN-=.33?»M*PP**(1  ./3.)*KFST1**.5/31 

- CALCULATION  OF  PC,  V  AND  H  AT  TPr.  kM<  OF  TrtF 

IPaNMTION  REGION 

.  5/1=4. *S 

PC  =  TBKS?1  ,  F’ S 1 ' A T  ,  N M S ) 

V  =  vc*S0kT ( 1  .-PC) 

»F.ST?l=PFFTu*6?l*V/Vr 

52  =52 1  +  2 ,  *r,s 
WO'  OJffTfcllK 

PC  =  TPP(S2,Pi>DAl  ,KPS) 

V=V0*SQKT ( 1 .-PC) 

Ke.S  1  2  =  PF.FTO*S2*V/VO 

JF(AHS(REST2-1.2Fn)  .  LF .  100.)  GC  TO  *0 

S2 1  =S21+  (S2-S21  )  /( PES'l  2-P  E-ST  21 )  *  ( 1 . 2E6-PFST2 1  ) 

S  2  t  =  S  2 

PE.5T2  1  *F<E'ST2 
52*521 
CU  TO  HO 
90  CONTINUE 

MT'JKM  =  0 .02  9o*M  *PK**(  1 .  /3. )  tRE.ST  ?**.H/S2 
N  i  l  =  2 

95  CuNiTMIfc; 


—  - - -- —  lUTfe.WP'Ol.ATTON  FUP  H  BETWEEN  THF  VAlUF.->  hT  TnE 

BEGINNING  OF  THE.  TPAtiS  IT  I  ON  RE.GIPiv  AN  U  AT  TnF 
E.rtU  OF  THE  TRANSITION  REGION 

HsHLAM+(HTHHR-mLAm')/1  .F  6*(PE  5-?.F5) 

WRITE! I W  ,  4  4  0  ) 

GO  TO  .105 

loo  continue; 


- TURhULF  NT'  FLOW  KLGIUN 

H=.0296*Kl*pP**d  ./3.  )*RFS*,*.B/S 
WRITFCIW,450) 

105  CONTINUE. 


- - - - CALCULATE  LOCAL  CONDITIONS 

Tb*R*TA*d.-'Ph**.b)*f  1. -PL/P A  J/CP/J 
CPB=. 2365+7. 6* Th/10.**6 

EFF=H*(TS-TA-T?+Tb)*SH/(WH*CPB*(Tb-TS) ) 

- — - ---OUTPUT 

5IN  =  S* 12. 

HPHsH*  3600 . 

TSF=TS-459.6R8 

THE*Tb-459.6Pb 

WKlTFdw  ,500)  SIN,HPH,T.SF,Tt?F,FFF 
DQ=H*DS*dS-TA-T2  +  T5) 

DTP=DG/(WB*CPb) 

TB=TB-DTB 

F:FFiC(I)=KFf  +EFFTCd) 

SINCH(I)=SIN 
S=S+DS 
150  CONTINUE 

160  Continue; 

WRITKCIW, 1300) 

DO  165  1*1, N 
EFFAVG*fcFFTC(D/NCASF 
WRITE l Ik , 1 1 00 )  STNCH(I),  EFFAVG 
165  continue; 

WPITFdW,  1200)  NCASE 
170  CONTINUE 
200  CONTINUE 
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300  »■-()(. v a 't  noe  p.?) 

400  F  ijP  V  AT  {  40  A  2  ) 

410  cuPV61 ( 1 H 1 ,40m?) 

420  FOR^Al  (  bX  ,  '  **  **•) 

430  FnRP  AT  C  5X ,  '  **  LAFjUAW  **') 

4 40  h  u»*.AT  (bX  ,  '  **  THAUSlTirifi  *<") 

4bO  FOWfA.T  (  5X  '  **  TuPbULFMT  **’) 

SOO  FfjBMAl  (4F15.?,F15. 3//) 

600  F  0R«A7  ( //2?X  ,  >  LtAOTN'G  FDOF.'  RADIUS'  ,  BX ,FS. 3 ,  1  FT'/ 

122X  ,'HEATtD  ARFA  F-FP  FOOT  SF A N '  F 7  ,  3  ,  5  b.OF'f'/ 

2?  2  X  '  1, 1 P  HlOHLIGrtT  DIAMETER'  ,4X,f6.3,  1  Fl'//) 

700  FOBVAXC /22X  ,  '  Al/miliiF  '  ,  15X,F  7,  }  '  F+'/?2X  ,'AVb  l  F  >  P  1  ,  1  7  X  ,  F  b  .  1  ,  1  F  1 

1/22X  ,  '  AIPSPr.FO  1  ,  1  7  X  ,Fb.  1  ,  '  KTAS'/2?X,  '  Ft  L  F  F  C  Fl,n*  PFH  r  G  1  ,  d  X  , 

2F4. 1  '  L  F  /  F  1  F-"  /  2  2  A  ,  '  P  L  F  t  D  Tr.MP  '  .  1  5X  ,  F  b  .  1  ,  '  F  '7  ) 

BOO  FUBv.A]  (///10X  ,  'b(I'l)  '  ,7X,  '  (P  FFh  hB  )  '  ,  3x  ,  '  (6F  I  N  T  f  R  ,  F  )  '  ,  /  X  ,  '  C  >- 1  F  *■ 

ID  TF.HF  ,F)  '  ,7X,  '  tFF  '/) 

900  FuBMA'Hlw  ,'S=',  Elb.B.SX, 'PC=' ,F.i*>.B,SX,  'TFiFIA  =  ',Fb,2) 

1  100  F.dPpAl  (Fife. 2, F  lb. 31 

1200  FuRnA'J  (1Hu,10X,  'Hu;.-bFH  OF  CASF.S  AVFFtAi,Fl>  s  ',13) 

1  300  FuP'A'l  (lHl.lOx,  '  M  T  0  1  '  ,  HX  ,  'FFF(AVG)  '//) 

1  400  FuPMAi  (  16FS  .  2  ) 

C A T.  L  F.X1T 
F  f.C 


A. 4  Function  TRF  Listing 


13o 

140 

150 


FUNCTION  1 Wk(  pSuAT  ,  I  MAX  ) 
riWtN.STGM  us[iA7Mb,2) 

Hu  i4u  K=  1  ,  7 

if  r s*i  ?.-pAuai  rK  +  i ,  n  n -  ,  no,  no 

Ci) '■■•I  TMlf: 

H’f  p5HA  1  (K  +  ^  ,  n >SOAT(  m  ,  I  )  }  140,  14'C,  IbO 

f  fjf’T  I  i-’i'F 

rOMTwi'fc. 


PC  =  F5DAT(F  ,  7  )  4  (**l?.-PftOAl  CK , 1) ) * (  f  PSUAl  (  K4  1  ,  /I) -Par  ul  (  K  ,  ■?)  1  / 

1  (PSLiA  J  (H  +  ]  .  1  r  (  K  ,  J  )  1  +  fS*l  2.-PSDAl(m  ,  1)  j/lPM'Al  IM?  ,  I  )  - 

2F.sr  AT(K,  1  )  )  *  (  ( PSPrtT ( fc ♦? , 2 )-pSL A i  (F+l,V))/(.Pf*CAT(K  +  ?,l  )  « 

3  PM' A  T ( K+ 1 , 1  ) 1-f  KSHATf K+!  , 2)-PSPAT (*,?) )/ (PSOAT ( pvt- 1 , l )-PM)mT(K 

4  1 


TPF  =  PC 
p  fc  tup  I-. 
FUf: 


1)1) 


A. 5  Sample  Input,  flata 


ill*  0  £5S;  rT*?^5  L7P  liEhl‘  AT>  ChANnM. 

V-  .v>  .:i3-!;!n-i;3D  3-°,  3-s  4.o  4.5  s.0 

3lb4  I  3V1, ^  lf’,J?AT  i  ft*-  j  oZ  uPbFp  [jp  7  *  *-r  K  *7  r  2“  1  '  ,'’f;’ 

ft S.4* *3f.o  ,  t'V.  tr^f:rm' 

184.  J*?.  264  ?;a  fl'  ?;*  3.0  3,5  4.5 

?1A3  ,  cni-Fj r-nPAT J J iipJimLic4*  9Ai 

To*  0  I*9,  ,  •?**.!  12  U7f-PF55JL-  f  m  c^ofu 

.?»:■  J  }?  ilif  ,*«;  ,„V  *«S’*1,5  .  s 

,  CGNFiruiUTlf.;*,  r?2* 

oil-  oj-2,  a  !<!•"  ‘4ft !<> 

20h.  21b.  ?32.  2bJ.  jh'  j?aS  }  h°  J  r  *  4.5 

17Sb  -  ^Pl^UBATjniJ  jit'ulXrZ  /?£•  »0fi. 

o,o- ».;c-?5„  j«?-6-  <w  ”*  tc“tr»'> 

15§.  1*2.  175.  ifts.  i^U  £*5  3.0  2.5  4.5 
1740  7  C0NFjCuA.4Tir,il  hZ'upIU-  ,  {?*  74,* 

0.0  0.5  J.o  j  2**  (  *  V?  344.$,  CLhbhClH- 

172‘  180.  143.  ?0U.  if,?.';  2;5  3,5  4,5 

1/4R  ^  CnwHGIit-ATiOiv  1  5/i*,/pL>jri  ,  ff*  ^8* 

1  •  0.;h-?  c  i’?.’  i'.s«  «*■ 

'*5-  i»«.  703%,'C  ,0"  ,hS  |i«  3,5  ,.b 

1738.  -h.^0tF,w5*5',n"  Jm'iiw-c-;  iff'  ,T. 

0.0  o#5  1  0  4 2*2.  39Q.7  or.vFf.riff, 

19§.  2fl«.  224?  237.  1  ^U  ,?$*  ,  h°  h 5  4*^ 

1740  •  7  C'DIvF  l  GUP  AT  If',’:  ]U%iPPFC  j  *  *05. 

2  •  o.;7-,  „  1*2- 7,  £  »**  cc^tnri. 

2«-  ,14.  73n.  jJj?  ,|5“  ,J-5  ,3.<.  5.5  4.5 

1740  7  CPNFJdJkATl  ,*  9  A  *UPPFP  fc*  J  <'*• 

1  8-  O.s'.'f  0  l4!-4  I5,S«  457^f  “f  Ce*»»C1lfl- 

775.  734.  ?J}^  ,}*«  2-5  i.O  3  5  4.5 

279a  1  9*’  UpCfp  }rl'  12 1  ,* 

0.0*  0.5  2i60  j*s,2->  350?8  7Af‘  CG88tClfU 

136.  139.  156?  ,Js;  ,/;V  2  3  0  3  5  4.5 

281ft  .  CONFH.jRATinw  •>**'  Ut>»t  ,  S^‘ 

S*  o.^°i  a  }^>.  0*VS  "Ko';7*  ^mCTFO 

152.,,  156.  ,174.  ,|^?  3.0  3  5  4.5 

2294.  -20?§NFIf9ftA3iPA'r?}fQUP^  COftfttrm 

0.0  0  5  J  /»  J’D*’  1  j£  ,  0  9  401  1  r  tUKKtCff.O 

158?  »ir?  ifaf  ,k?  ?55  j  o '■’•]  5  4 .5 

0.0*  0  i20?Bri  }^*  ii%  P  lb iL¥  T4f  cnKAtrm) 

i»«?  1S5?  2^  2 U:  ^5°  ?5S  h°  l>5  4.5 

2775  _CnNFX(;URATinw  23AUPPH  ,  ru’  y?/,. 

.  !sf->,  17 ftUP^,,;{1'  JS>  cc8»fClflJ 


• 

*  ^  ^  f  C  ft  F  C.c  J  f  f, 


4.5 
74. 

'J4f  LCPhtCl  H 

4.5 
68. 

rrppf.c/’f  u 

4.5 
92 . 

7Af  CCPbtClTU 


1AF  CG^OFC'lM 

4.5 

121. 

7  4f,  CCFftfCtfeO 

4.5 

b6. 

T4b  CCPPtCTF  0 

4.5 

ft8. 

l'4f  CCfcfrfc.CTf.Xi 
1‘4f  CCRPtCCTKo 


0.0*  0,5  i  0  P?*1-.  12.03  451,1 

180.  185.  20fi?  226?  16BU  Iv  |l°  |.5  4.5 

274U  -20?5f,K1?y4ATIn\24^*Upp^  L IP* 

20 5?  205?  ai«?  Al*\h«  Al*  i?o5,*37.s  4.5 

.  9ft 


t.  k  f  i  c  i  k  >  r  y 


iAf  COfcfcf  CTfc.O 


A . 6  Sample  Output  Data 


FUuF'L.  36  FuGIM  IM.Ft  UP.  b'l.tFl)  Aik  CHANf-FL  t,f'r  h.I:.?  Cl 


LkrtDJf'.G  FLGF  HAT  Hus  '  0 .  J  2  5  t  r 

HKATkp.  A  r  I  A  kkk  FOOT  SPAN  0.55R  S.jFi 
OTP  HIGHLIGHT  U  I  A  h  t,  T  £  R  2.237  A  I 


CONF  IGUhA'i  UJH  1  (»A  UpPrk  LIP 

aM  1  f.ll  a 
hf'k  THMP  . 

A  1  h  GpFhl' 

HI  t,  FI)  P'Jill'-i  PPP  F  IjG  ]  f>‘ k. 
E>Lc.F,i)  1  k  * •"  P 


TaF  CfiRRFCTtl) 

3164.0  Ft 
1.3  F 

2  2  4.8  FT  AS 
j  4 . 9  l.k/i*T ( 4 
402.4  F 


MIN) 

S~  0.0000G000t+00 
**  CYLINDER  ** 

0 .  no 


S=  0.4lboo999F_-01 
**  CYLINDER  ** 
0.50 


S=  0. 833339998.-91 

**  TRANSITION  ** 

1.00 


S=  0. 12500100K  +  00 

**  TRANSITION  ** 
1.50 


S=  0.1666b800t+00 

*♦  TRANSITION  *# 

2.00 


(H  HK»  Hh)  (SKJN  TtMp,8  )  (PLPr.n  T,,:4*,n 
PC=  0. 1  onoOOOOf +01  THf.’I  A  =  O.j" 

4  3.0  0  1*4.00  <♦  0  2  .  4  <  • 

PC=  0.9b949o57F+00  J'HRTA  =  19.10 

42.59  1 92.00  392.54 

RC=  A  .  1  299bl.>  30F  +  0g  THETa  =  3k. 20 

21.74  204. Ou  382.32 

PC  =  -0. 14700172E+01  THF.IA  =  57.30 

39.10  223.00  370.70 

PCs-O. 12 999b 4 8F +01  THKTA  s  76.40 

47.34  160.00  3  6  h . h  4 


S=  0.208  3  3500  8. +  00 

*♦  TRANSITION  *♦ 
2.50 


PC s-0, 93949 0  22E  +  00  ,  THfc.7  A  =  9b. 5O 

52.81  100.00  355. h3 


Ss  0. 25000200 H.+ 00 

**  TRANSITION  ** 
3.00 


PC  =  -0. 8099961  3E.  +  0U  H'FTA  =  11a. faO 

59.56  94.00.  349.41 
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!■'  r  K 

0  .  c<>  t 


0  .  C  1 


0  .  a  2  2 


0  .  *A  *  I 


u  .  e  2  r. 


0  .  i  3  p> 


t  .  36  j 


SCI'.) 


f  F  F  f  a  y  r,  i 


o.OO 
n.SO 
1  .  C 
J  .50 

?.oo 
>  o 


*  .  00 


**  u  “  y  f  u, 


0.^7 
0  .  *>  P 
\44', 
0.  »S7 
0  .  ►  0<5 
0  .  1  1  7 
'  .  3  ?  n 


*■  OA  S‘S  A  V 1  i»  J o. F D 
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APPENDIX  B  -  ICING  ANALYSIS  PROGRAM 


Flow  Chart 


Flow  Chart(Cont-) 


END 


B-8 


Input  and  Operation 


B.2 


MODE  -  1 

This  mode  of  operation  reduces  experimental  data  to  check  for  accuracy. 

An  expanded  output  shows  results  of  intermediate  calculations  in  the  point 
by-point  analysis. 


INPUT  PREPARATION  FOR  MODE  -  I 
CARD  1  ( 20A4 )  Title 

This  card  can  contain  up  to  80  columns  of 
alphanumeric  information.  It  will  be  used 
as  the  heading  on  the  putpiit  data  sheet 
(20A4  format). 

CARD  2  (1F1Q)  Mode 


Enter  a  "1."  in  column  1  for  Mode  I 
operation  (F10. 3  format). 

CARD  3  (8F10)  Nacelle  Geometry. 

This  card  has  8  fields  of  10  columns  each 
to  contain  the  variables  C,  TOC,  RLE,  DIA, 
SH,  SW,  COO.  (F10.3  format). 


Col 

umn 

1 

C 

Col 

umn 

11 

TOC 

Col 

umn 

21 

RLE 

Col 

umn 

31 

OIA 

Col 

umn 

41 

SH 

Col 

umn 

51 

SW 

Col 

umn 

61 

CDO 

Col 

umn 

71 

OCD 

Chordwise  length  of  nacelle  in  feet. 
Thickness  ratio  (maximum  nacelle  thickness/ 
chordwise  length)  non-dimensional. 

Radius  of,  curvature  for  the  leading  edge  of 
nacelle  Inlet  lip  In  feet. 

Nacelle  Inlet  highlight  diameter  In  feet. 
Streamwlse  length  of  the  heated  surface  In 
feet.  (Heated  area  per  foot  circumference). 
Aircraft  wing  area  In  square  feet. 

Drag  polar  term  (C^  or  CD,„ter(.ept)_ 

Drag  polar  term  (k  or  Cn  ) 

°slope 

(Cq  ■  CDO  +  (DCD)C^  Drag  polar). 


CARD  4  (15F5)  Pressure  Coefficient  Table. 

This  card  has  15  fields  of  5  columns  each  to 
contain  streamwise  distance  measurements, 
from  stagnation  point  to  point  under  analysis, 
in  inches.  The  values  shall  begin  at  zero 
(stagnation  point)  and  increase  at  any  de¬ 
sired  interval  until  the  aftmost  heated 
region  is  reached  (F5.0  format). 

CARD  5  (15F5)  Pressure  Coefficient  Table. 

This  card  has  15  fields  of  5  columns  each  to 
contain  pressure  coefficient  values  for  the 
corresponding  points  on  card  4  (F5.0  format). 


CARD  6  (10F7)  Efficiency  Profile  Table. 

This  card  has  10  fields  of  7  columns  each 
containing  the  streamwise  distance  measure¬ 
ments  from  the  stagnation  point  to  the  point 
under  analysis,  in  inches.  The  values 
Increase  at  any  desired  interval  until  the 
aftmost  heated  region  is  reached  (F7. 3  format) 


CARD  7  (10F7)  Efficiency  Profile  Table. 

This  card  has  10  fields  of  7  columns  each 
containing  the  channel  efficiency  values  for 
the  corresponding  points  on  Card  6.  (F7.3 

format). 


CARD  8  (20A4)  Test  Title. 

Test  number  and  any  pertinent  information  up 
to  80  columns  (20A4  format). 
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B.2  Input  and  Operation  (Cont.) 

CARD  9  (3F10)  Flight  Parameters.  This  card  has  3  fields  of 

10  columns  each  to  contain  the  variables  ALT, 
W,  VKTAS  (F10.3  format). 


Column  1  ALT 
Column  11  W 
Column  21  VKTAS 


Altitude  (MSL)  of  aircraft  in  feet. 
Weight  of  aircraft  in  pounds. 

True  airspeed  of  aircraft  in  knots. 


CARD  10  (4F10) 


Column  1  TA 
Column  11  DDM 
Column  21  LWC 

Column  31  HORIZ 


Icing  Conditions. 

This  card  has  4  fields  of  10  columns  each  to 
contain  the  variables  TA,  DDM,  LWC,  &  HORIZ 
(F10.3  format). 

Ambient  Temperature  (in°F). 

Water  droplet  median  diameter  in  microns. 
Liquid  water  content  of  surrounding  air  in 
grams/ cubic  meter. 

Horizontal  cloud  extent  in  statute  miles. 


CARD  11  (3F10) 


Column  1  WBM 
Column  11  TBIN 
Column  21  NIP 


Bleed  Air  Data. 

This  card  has  3  fields  of  10  columns  each  to 
contain  the  variables  WBM,  TBIN,  NIP  (F10.3) 
format) . 

Bleed  air  mass  flow  rate  per  engine  in  lb/min. 
Bleed  air  temperature  in  °F. 

Percent  engine  RPM. 


In  this  mode  a  dummy  card  with  "33."  punched  starting  in  column  1  should 
follow  card  11.  This  will  loop  the  routine  back  to  read  In  a  new  set  of 
test  data.  Cards  8,  9,  10  and  11  followed  by  a  "33."  card  should  be  re¬ 
peated  for  as  many  test  cases  desired  up  to  100  tests. 

One  blank  card  following  a  dummy  "33."  card  will  allow  the  Input  of  new 
mode,  nacelle  geometry  and  pressure  and  efficiency  data  cards.  Two  blank 
cards  following  a  dummy  "33."  card  will  terminate  the  program. 


B.2  Input  and  Operation  (Cont.) 


Mode  II 

This  mode  of  operation  is  for  icing  prediction  under  given  flight  and 
icing  conditions.  Intermediate  information  or  calculations  are  not  shown. 
INPUT  PREPARATION  FOR  MODE  II 


CARD  1  (20A4) 

CARD  2  (1F10) 

CARDS  3-7 

CARD  8  . 


CARD  9  (4F10) 


CARD  10  (3F10) 


Title. 

Prepared  same  as  for  Mode  I  operation. 

Mode. 

Enter  a  2.  in  column  1  for  Mode  II  operation 
(F10. 3  format). 

Same  input  as  described  under  Mode  I 
operation. 

Flight  Parameters. 

This  card  ha',  3  fields  of  10  columns  each  to 
contain  the  variables  ALT,  W,  VKTAS,  (F10.3 
format).  These  variables  are  input  the  same 
as  for  Card  9  under  Mode  I  operation. 

Icing  Conditions. 

This  card  has  4  fields  of  10  columns  each  to 
contain  the  variables  TA,  DDM,  LWC,  &  HORIZ 
(F10.3  format).  This  data  is  input  the  same 
as  for  Card  10  under  Mode  I  operation. 

Bleed  Air  Data. 

This  card  has  3  fields  of  10  columns  each  to 
contain  the  variable:.  WBM,  TBIN  and  NIP 
(F10.3  format).  This  data  is  Inf  the  same 
as  for  Cc.'d  11  under  Mode  I  operation. 


In  this  mode  up  to  100  sets  of  cards  9  and  10  may  follow  a  Flight  Parameter 
card.  A  dummy  "33."  card  following  a  card  9  and  10  set  will  allow  input  of 
a  new  Flight  Parameter  Card  8.  One  or  two  blank  cards  following  a  dummy 
"33."  card  has  the  same  effect  on  the  program  as  during  Mode  I  operation. 


non  non  non  non  non  onnn 


P.3  Program  Listing 


PROGRAM  ICFOFT  . ...»  .6/22/79 . . 

ENGINE  INLET  LIP  ANTI-ICE  ANALYSIS  _ 

DETEPMIuATIUN  OF  ANTI-ICING  PERFORMANCE 


REAL 
PEAL, 
PIvFNCTON 
01 Mt NS  I UN 
01  PENSION 

dimension 


J  ,  K.KO,  K1  ,  K2  ,'L  ,  L0LS,LWC  ,  MU,MW,M*AV,M*tMP  ,  MWLE,  .*1  K,  IMFPH 
MODF 


P  S  0  A 1 (lo, 2) 

EFPAT( 1 o, 25 
T  1TLE  (  2(3 )  , 

TEST (20) 

OPEN (UNIT  =  2 , NaME= 1  TCI  OFF. DAT ' , TYPE=’OLD' ) 
0PFN(UNIT=3,NAME=' ICFOFF.uUT' , DISP= ' PR  TNT ' ) 


20 


25 


•CONSTANTS 


I  H  =  2 
Tw  =  3 
NPSs 1 6 
NESslO 
G=32. 174 
Pl=3.1416 
J=77P. 

P  =  5  3 . 35 
DS=. 0416667 
DTSJN=20. 


•Inputs 


1  =  1 ,NPKT , 
1=1 ,NES) , 


PbAD ( IP , 300 )  TITLE 
DO  277  J J= 1,20 
PEADfIP,310)  MODE 

READdR,  310)  C,TOC,HLE,DIA,SH,SW  ,  CDCl ,  DCb 
IF(hLE  .LE.  0.0)  GO  TO  2P0 
PEAD( IP , 1 100)  ( ( PSD AT ( I , N ) , 

PEADf IR, 1200)  ( ( FE  DAT ( T  ,  N  )  , 

DO  276  1 1  s  1 ,100 
X  M  =  0 . 0 

I F ( MODE  .NE.  1,0)  GO  TO  IS 
READ ( I R , 300 )  TEST 
15  CONTINUE 

READ  dR  ,3 1  0  )  ALT , W , VKTAS 
IF ( VKTAS  .LE.  0.0)  CALL  EXIT 
IFCVKTAS  .GE.  1.)  GO  TO  20 
XM* VKTAS 
CONTINUE 

WRITE(Tw,330)  TITLE 
IF (MODE  .NE.  1,0)  GO  TO  25 
WRITE ( I W , 320)  TEST 
CONTINUE 


N=t ,2) 

»--l,2) 


■CONVERSION  OF  INPUT  DATA  TO  INCHES 


CIM=12.*C 
RLE 1H=12. *PLE 
DIAIN=12.*DIA 
SHIN=12.*SH 

WRITE  dW,  400)  C1N,T0C,RLEIN,PIAIN,SHIN>5W»CD0,DCD 
WRITE ( IW , 500 )  W 
WHITE ( IW , 600 ) 

DU  275  IJsl.100 

READ(IR, 310)  TAF,DDM,LWC,HORIZ 
IF (TAF  .GT.  32.)  GO  TO  276 


•CONVERSION  OF  UNITS 


SL»SH/2. 

N=SL/DS+ 1 
TAsTAF*459.68B 
DD=DDM*3. 281/10. 4*6 


•CALCULATION  OF  BASIC  PARAMETERS 

B-13 


nonooinaon 


.  As5  <HT ( 1 .4*G*R*TA) 

.LF.  0.0)  go  rc  40 
VKtAS=X^*A/i.b»7H 
40  V  f V  *  T  A  S  *  1 , o  A 7  8 

rp=. 23*5*7. 6*TA/10.**fc  „ 

*»0s7.475*1A**1.5/(TA*216.)/10.**7 
t‘l  =  (  4.722+.ub^4,*lA)/l0.**7 
P  R  r  C  P  *  MU /  K 1 

PA  =  2  11*  .2  1  /EXRf  AL1/  ( 277 1 0 . - . 09877 4* ALT  J  ) 

RpC=Fa/TA/P 

PtFTf".spKr*VC»/**«l 

PtI  r  =  RFFT(j*2.*RLF 

X  M  r  V  0  /  A 


•CAt  CUlATI ON  Or  BLEED" A  T W  DATA 


DRrpMO/2 .  *Vf>*VC'/G, 

CL=» /DP/S* 

Fi<=CCDfJ*DCr>*CL*CL)*DP*S»/2. 


•iNStRT  calf  FOR  ENGINE  Hr RFURR ANCfc 
PkOGP  a K  it,  uPT  A  J  N  BLEED  A1P  FLO* 
KATfc  AfiL)  TFMPEHATUPE  C)»  PEAb  THIS 
DATA  from  input 


PEADf IP, 310)  «BK, IMF, NIP 
T  6  T  \'  =  1  B  r  ♦  4  5  9  ,  o  8  8 
*ri=*e»i/60./Pi/niA 


•CALCULATIuN  of  impingement  data 


IFfLWC  .GT.  0.0)  GO  TO  60 
P«<Lc  *0. 

GO  10  90 
60  CONTINUE 

PEC=BEFTC*DD 

IF (RED  .GT.  200.)  GO  TO  70 
68  LOIS*. 98-. 1 J4*AL0Gf RED) 
go  ro  71 

70  LOl S=.74-,OH87*ALUG(PeD) 

71  «  2s  .  1  08*VO/Ml1/C*DD*Dn*G 
FO=K?*LOLS 

IFfKO  .GT.  .01)  GO  TO  80 
IFfKO  .LT.  .004)  GO  TO  77 
75  F**.0873*( A LOG fKO)*5.522) 

GO  TO  85 
77  EM*0. 

GO  in  Kb 

80  IFfKO  .GT.  .*)  GO  TO  83  ~T“~T_ 

EMa.08T,31»(2.>.4J42*ALOG(KO))**l .55 
GO  TO  85 

83  WRITE ( I  * , 900 )  KO 
GU  TO  275 

85  WNeVO*LRC*TUC*C*EN*.623/10,*M 
RrfAVswH/SH 
M*LE*2.*MNAV 

90  continue 


•POlwT-PY-POINT  ANALYSIS  OF  HEAT/MASS  • 

TRANSFER  RELATIONSHIPS  STARTING  STAGNATION 
PuInT  AND  INOING  AT  AFTMOST  POINT  OF  HFmTEO 

REGION 

HEATED  REGION 


•INITIALISE  CONDITIONS  AND  START 


SapSDAT (l»1)/l2. 

NN*1 

tb*tBin  _  -  n 
ThFI N*TB IN-459, 6R 8 
DRR*0.0 


nnno  noon  ooo  non  non  non  nno  mo 


TOT  1 M  =  0 . 0 
DO  2  53  1*1  tN 

IFCS  .It.  SL)  GO  TO  ICO 
S  =  SL 

100  SP  =  S/SL 

- CALCULATION  OF  LOCAL  HEAT  TRANSFER  C'ltFr'ICIfcM 

IF f I  .01 .  1 )  GO  Tu  1 15 

- - CALCULATION  AT  THE  STAGNATION  POINT 

PC  =  PSi)  AT  (1,2) 

PL=PA*DP 

THFTA=S/PLFa57 . JO 
PKSsO.O 
CO  TO  120 
1  1.5  CONTINUE 


......... - - - CALULATION  am a y  from  the  stagnation  point 

PCslPP(S,PSi)AT,NPi>.IM) 

.  THFI  J!=S/.PLE*57.  JO 

V=V0*50KT( 1 .-PC) 

PL=PA^PC*DP 

PtS=REFTO*S*V/VO 
1 20  CONITnUE 

IF  (THFTA  ,GT.  25.)  GO  TO  125 
- . - - - CYLINDER  REGION 

H=. 57 *K1*PR**.4*P£LE**.S/PLF* Cl .-(THETA/90. )**3) 

GO  TO  1 H  5 
125  CONTINUE 

1  F  (  RFS  .  GT.  2 .  £5 )  G‘J  TO  135 
...... - - - - INmW  FlOM  REGION 

Hs. 332**1 *PR**( 1 ./3. )*RES**.5/S 

GO  TO  IBS 
135  CONTINUE 

IFCPFS  .GT.  1 . 2t6 }  GO  TO  180 
- --TRANSITION  FLOP  REGION 


S 1 1  =  S 

PfcSTU*KfcS 
S1*S1 1 *DS 
GO  TO  (142,17b) 
142  CONTINUE 


NN 


CALCULATION  OF  PC,  V  AND  H  AT  THE  BEGINNING  OF 
THE  TRANSITION  REGION 


PC*TRP (51 .PSDAT.NPS) 
Y«VO*SOPTil .-PC) 
REST!*REFTO*Sl*V/VO 


REST! *REFtO*S l * 
IF(ABS(RESTl-2. 
511*51 U(S1-S«1 
SllsSi 

REST11*REST1 
S1*S1I 
GO  TO  142 
CONTINUE 


»V/VO 

i OEb )  ,LF.  100.) 
I  )  /  (RESTl-REST 1 1 


GO  TU  150 
)*(2.0Fb-RESTll ) 


HLAN*.332*K1*PR**( 1./3. )*RE5T1R*.5/S1 


CALCULATION 'OF  PC,  V  AND  M  AT  THE  END  uF  THE 
TRANSITION  REGION 


S21 *4. *5 

PC*TRP(S21 ,PSDAT,NPS) 
V»VO*SQRT(i.-PC) 

8-15 


onnno  nnonn  nnr>  ono  nnn  noaon 


PEST2l=REFTU*S2l*V /VU 
S2=S?1+2.*DS 
• lhO  CONTINUE 

PC=TRF  fS2,PiiDAT ,’iPS) 

V=vu»SORT( 1 .-^C5 
WE51  ?s:PEFT0*S2*V  /  <jn 

IF(ABi>(RFST<-i.2F8)  .  IjE  .  100.)  GC  TO  HO 
S2I=S2l4(S?-S21  )/(KESi?-RES[2n*n.2£6-RFST21) 
S  2 1  *  S  2 

PEST21SRFST2 
S2=S2I 
GO  m  160 
170  CONTINUE 

HlURB  =  0.029t>*M*PR**(  1  . /3.  )*PES12**.8/Si 

jy  ^  —  2 

175  CONTINUE 


•TnTERrOl,A  1  TON  FOP  H  R£TrfKfc.N  Tnfc  VALUE*  AT  T-.E 
»tG!SNT,-G  OF  THE  TRANSITION  REGION  ANO  AT  Trt^ 
ri'.n  OF  THE  rPANSITTO'J  REGION 


H=HLAM>(H1UKB-HLAM /I .K6*(WES-2.E57 
GO  TO  1 85 
180  CONTINUE 


- - - TL'BodLKrtT  FLC1*  REGION 

H-.0?96*K l*pR**C 1 . /3. )*KFS** ,8/S 

- - 0  ALCULATF  I.OCaL  I  NPINGEbFNT  (4HUnBACK)  PATE 


185  IF'SP  .GT.  .  3  3  3  3  )  G'l  70  1  HP 

**IMPrM*Lfc.*( 1 .-.3S5A(3.*SP)*»1.75) 
GU  TO  180 

1H8  PwI8P  =  P-Lb*l  .177*11  .-SPj**1  .  6 
189  PHr*»El*«iP  +  DPt» 

TUTI*  =  TOTI*-*»*IhO 


CALCULATE  T2.T5 


PLsPA+PC»OP 

T2  =  VO*VO/(2.*CP*0*.J)*(PH**,54CP*NW/H) 
TSe(l .-PR**.5)*R»TA/CP/J*(1 . -PL/PA) 


- -- - calculate  ts,ti,T3,t4 

CASE  I.  NU  EVAPORATION 

K*0. 

CPR*. 236547. 6*TB/10,**6 
EFFaTHP  (  S  .F.FOAT,  N£S) 

200  TS*  (  WB*CPB*£FK*TR4H*SH*  (TA*  (  1  .  4NN/H  )  ♦!  2-T5  )  )  /  C  NB*CP»*bFK4H*SH*  (  1 
lNP/H)) 

T1*(TS-TA)*(1 .4MH/H) 

OIN  **B*CPB*KFF*(1«-TS) 

Q0UT*H»SH*(T1 -T24T5) 

205  IE  (TS  .LE.  491. 6P8)  GO  TO  240 
IF(*W  .LE.  0.)  GO  TO  240 


210  TSMOl.bfi* 

DTSsniSIN 

212  T1*(TS-TA)*(1 .4NN/M) 

L*1 348,21-. 562*15 
ES»PP ( TS) 

TJ*0,b22*L*ES/CR/PL*144. 
FA*RP ( TA ) 

T4a0.622*L*EA/iP/PA*»44, 

OI*»*MCPB*EFF*lTb-TS) 


■CASE  2.  NITH  EVAPORATION 

KINO  KOUIL1HIRIUN  VALUE  uK  Ti» 
SUCH  T*AT  OIN*OOUT 
BY  ITERATING  BETWEEN  492  AnO  TH 


B*16 


onn  nnnn  non  non  on  no  non  non 


B . 3  Program  Listing  (Cont.) 


C 


OlJI'TsH*SH*m-T2  +  T3-T4  +  T5) 

IFCLIN  .Lb.  Ci OUT )  GO  TO  219 

IF  UJN.GT.OOUT,  INCREMENT  TS  BY  OTS 

TS=TS+DTS 

Tf  as  ,LT.  TQ)  GO  Tt)  212 

WWTrrci*.  700)  ALT, VKTAS.TAF, NIP.NBR.TBFIN, CD« ,L*C 
GO  TO  275  ' 

2)9  TFC1S  .GT.  491. 68b)  GO  TO  225 

If  UIM.LT.QOUT  BUT  TS.I.f  .491  ,t>8d  , 

FOoJLIBIRIUy  IS  »  491. 6»R 

FIND  EVAPORATION  RA1E  F  UR  Tz  =  491.  nbS 

720  7S=491.68» 

0INsWh*CPB*bFh  *(Tp-i'S) 

001*1  Sul  N 

Tl=(TS-TA)f (1 .4RN/H) 

T3eQOUT/f M*SH)-(Tl-T?-T4*T5) 

go  rn  230 

225  IFCGOUT/QIN  .lE.  1.001)  GO  TO  230 

IF  OIN.LT.OniJT  HUT  XS.GT.4  92, 

DECREMENT  TS  BY  SMALLER  AND  SMALl  r.9 
AMOUNTS  TCJ  OBTAIN  CON  VFRGhNC’F 

TS  =  TS*f)TS 
DTSsDlS/3. 

TSzTS+OTS 
GO  TO  212 

IF  OIN  s  OOUT,  CONFIRM  THAT  KVAPOr A  l  TLU 
,  RATfc  DOES  NOT  EXCEED  1MRI NGft.Nh.iVT  RA  TF 

230  Fr(T3“T4)*W/fD*NW) 

IF(K  .Lt.  1.)  GO  TO  240 

IF  EVAPORATION  RATE  EACr.FDS  l  HPI  NGEAFN' 
RATE,  FIND  NFR  TS  FOR  FULL  EVAPORATION 

T3=f 44L*M*/H 

TS=(  W»*CPh*EFF*TR  +  H*Sh*(TA*(  1  .♦MW/NJ+1  1 5  )  )  /  (  »t>  *CPH*EFF  +  H* 

1SH*(1.-»MH/H)1 
0IN5WB*CPb*FFr  *CTa-TS) 

TI*(TS-1A)*(1.4HW/H) 

00UT*H*SH*(Tl-T24T3-T44T5) 


....... - .....  CALCULATE  ELEMENTAL  RUhBACK  RATE 

240 


CALCULATE  DROP  IN  TB  AND  GO  TO  NEXT  PUTnT(SKGMFNT) 


DOeOOUT*DS 
DlPsDO/C Wf*CPH) 

T8»TR-DTB 
SI N*S* 1 2  » 

TSF»TS-459.6R& 

HPH*H* 3600 , 

DPR H*DRB* 3600. *PS 
DMWJHaN«IMP*3600.*DS 
TEC  MODE  .  N£.  1.0)  GO  TO  25? 

CONTINUE 

WRITE (I *,1000)  SIN,PC,EFF,TSF, DMW1H,K, DRBN, HPH 
CONTINUE 
S*S4DS 
CONTINUE 

.............  CALCULATE  TOTAL  IMPINGEMENT ,  EVAPORATION 

AND  RUNBACK  PATES 

RBPH*DRB*3600,*DS 
I«PPH*TOTI**3600.*DS 
cvaph*impph-rbph 

IFCINPPH  .GT.  0.0)  GO  TO  262 
260  PERC*Of 
GO  TO  265 

262  PERC*EVAPH/I*PPN*tOO. 

.................  CALCULATE  PUNPAC*  CROSS-SECTION AL  ARM 


250 

252 

253 
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non 


265  TIME=hnRIZ/VKTft6/l .151 

PHS  =  nRB*3600.*Do*llMF.  / 6 2. 4*  144. 


OUTPUT 


,  khrL^n*,'rt!r°(l  ALT  .  VKTAS,TAF,MP,*RM,TBFIN,DrM,LWC,HPHT£,IMPPb, 

X  fc  VArM  >  pK  «CfKr\PHfWP£> 

275  COMITMJt 

276  COMTNUfc 

277  CONTIMJr. 

280  COMM  vtUrJ 
3o0  Ffjpfc  A’j  (  20A4  ) 

310  FURBA1 (WF10.3) 

320  FuR*AT(lh0»b*,20A4//) 

330  FOP^A'i  ( 1 H 1  ,  19A.20A4//) 

400,o9S’?nT(^<4,^CfcVL?  LKNGTH'  r  F 1 «  .  2  e  *  IN*/23X.  •  THICK  NFSS-TU-LFNG7H 

1RAT10' .F8.3/23X, 'LEADING-LDGF  P AHl MS ' , F 1 3 . 2 , {  IN'/23X, 

i«kIP,AJX?£fcTt:R'  f  F?0 . 2  ,  •  IN  ' /23X  ,  'STREANwISE  SURFACE  •  /2  3x  ,  '  LENGTH 
A*  V.FH  1  Kl  •  /2  3X  »  '  W 1 NG  AREA'  F23.2,'  Su  FT'// 

«:„„423X'  tRAS  RL'LAP  ...  .CO*  '  ,  F6.4,  '  ♦  •  ,  F  6 , 4  ,  •*CL*CL'/) 

500  FuPf-ATf  23X. 'PfclGHi  •  .F26.2,  T  LB ' / / /  )  ' 

60°  FpPHA'i(3X,*A[J  VtL  AMB  PEPC  rtLfc.Fl)  BLEF  0  DROP  l*C  HOP  TMPJNG 
1  LV£?,-r-PE>9  F'Jf  bK  RUNHF ' /1 3X ,  'TKVP  Ml  FLO*  TEVP  PIA  (CM/  EAT 
i  ,fc;AP  PATE  SECT'/'  ( FT  1  (KTAS1  (F )  *  ,  6a  * ' (  FPW  )  (F) 

_„3  (*IC^  .Cm)  (»<•)  fPPh)  fPPh) ' ,bX , ' (PPH)  (Sul U )  •  ) 

TOOiFORKATU.  8.p,2F5.0, Fb.o.E  7.2,F6.0,F  5.0.F8.3.5X, 'NO  CUNVFkGENCF  ON 

?R^Mnf9'!:7xn{,,'?«04F6X2»2F5*0:f'6«3'3F6-2.):!>.O,F6.2,F6.3) 

iiooS/ip°J^SsjiF5-1’3X’  "* ,l 

1200  F'GPtf  AT  C 1  0F7 . 3  ) 

CALL  EXIT 
FNC 


nnan 


B.4  Function  PP  listing 


FUNCTION  PP(T) 

- FUNCJIuN  PP  CALCULATES  THE  PARTIAL  PRESSURE 

OF  SATURATED  AIR  If!  FOUNDS  PER  SQUARE  INCH. 

ZTslOOO.O/T 

IF fT  .GT.  491 . 6oR )  GO  Tu  ]Oo 
A=19. 598997 
Bs-10.4310'5 
Cs-0. 27550  73 
D=0. 039404  *93 
E  =  0.0 
GO  ID  300 
100  CONTINUE 

7FC)  . GT .  671. 6n«)  GD  TO  200 
A  =  1  3.435290 
P=-5. 09884^4 
Cs*l. 689617% 

0=0.17829154 
E-0.0 
GO  m  300 
200  CON1IMJR 

A=lb. 625544 
B  =  - 14.2131  Ob 
C= 7.5567694 
D=-4«01515*9 
F=0. 71697364 
300  CONTINUE 

XP=A+H*ZT*C*ZI**2*0*ZT**3+F*ZT**4 

PP=EXP(XP) 

RETURN 

END 


MODEL  35  ENGINE  INLET  LIP  ANTI-ICE  ANALYSIS  -  12  O'CLOCK  H’SllION 

1  *  7.814  .058  P.  125  2.237  0.558  253.  .0^25 

0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0 

1.  .97  .13-1.47-1.30  -.94  -.81  -.77  -.81  -.92-1.09 

.0  .5  1.0  1.5  2.0  2.5  3.0  3.5  4.5 

0*627  0.696  0.445  0.897  0.609  0.317  0.32R  0.35C  0.350 

TEST  8 1  IN'  1A  -  OPT  CORRECT  ED 
3195.  13000.  236.8 

16.33  0,  0.  0. 

348.9 


3195. 

16.33 

11.98 

33. 

TEST  RU.4 
3181  . 
15.86 
14.94 
33. 

TEST  RUN 
3189. 
16.08 
12.14 
33. 

TEST  RUN 
3 1  T9 . 
15.37 
15.0 
33. 

TEST  RUN 
3123. 
15.80 
12.12 
33. 

TEST  RUN 

3150, 

16.44 

15.05 

33. 

TEST  RUN 

3135. 

1.48 

12.14 

33. 

TEST  HUN 

3160. 

0.71 

15:06 

TEST  HUN 
3160. 
0.50 
12.19 


2A  -  DRY 
1  3000. 

3bS .  6 

3 A  -  DRY 
1 3000. 

4o6.9 

4 A  -  DRY 
1 3000. 

,  0. 

399.3 

5  A  -  DR  Y 
1 3000. 

44$. 1 

6A  -  DRY 
13000. 

°i 

451.3 

7 A  -  DRY 

1  3000. 

0. 

352.2 

HA  -  DRY 
13000. 

0  , 

352.2 

■  9 A  -  DRY 

13000. 

0. 

399.3 


COHRECTE  0 
235.9 

0.  0. 


COKRECfEL 

236.8 

0.  0. 


CORRECTED 

238.4 

0.  0. 


CORHECTEL- 

235.4 

0.  0. 


CORRECTED 

234.2 

0.  0. 


CORRECTED 

231.7 

0.  0. 


CnHRECTED 

237.5 

0.  0. 

COKRtCIE  D 
232.7 

0.  0. 


TEST  HUN  10A  -  DRY  COHRECTFD 
3154.  13000.  224.8 

1.33  0.  0.  0. 

J4.95  402.4 

TEST  HUN  1 1 A  -  DRY  CORRECTED 
3193.  13000.  726.1 

1.91  0,  0.  0. 

i$.07  4S1.1 

33, 

fEST  RUN  12*  -  DRY  CORRECTED 
3187.  13000.  226.0 

1,23  0, 

14.97  453.2 

3f E3T  RUN  13*  -  DRY  CORRECTED 
1756.  13000,  142.6 

-8.95  0,. 


8,5  ^aseje  Input.  Data  {fw  \ 


test 

1740. 

-7.J4 

4.94 

*est 
1748. 
-6 .51 
12.09 

3fEST 

1738. 

-6.81 

15.46 

»fcfT 

-6.98 

12.01 

fesr 

i?;s* 

IS:04 

--Jest  i 

2794. 

-22.  h 

H:4* 

2.I1?1  ' 
-20. §8 

1!:‘* 

II-5’4 

- 

-20.16 

M:1* 

27?!”  » 

4i;*f 

720.il 

55: ’6 


5T  RUN  14A  -  DRY  ruR«?mn 
i  1  3000.  142.1  CTtD 

*  349.0 

'T  HUN  15A  -  DRY  CURRJ.'f  ri.n 
13000.  142?7  ftD 

0  . 

402. 

I  «•«  }f»  ;  dry  CURkfCTED 

0.  * 

399.7 

T  “"’'nioo.'  0"v1J5^Ecif  D 

450. 

r  HUN  1«A  -  Day  rrjRfai.  rru  n 
13000.  '142.4  CTtD 

452.6 

’  RU“l3000.'  DRl,^05HFC,tC’ 
35o.e 

RUHuooo.*  DRY,$HfcTtD 
3sS.j 

RUN13?0*.'  "‘-{^RRECTtr 

401.1 

R"\jOOO.*  “"If^ORRECTED 
4«1.3 

RU"l30oS.'  DBf,«»'>ECTiD 

461.1 

RU"i3ooo.*  I>R\S?R?FCIt0 

4$1 .7 


B.5  Sample  Input  Data  (Cont.) 


1*1. tT  LIP  ANTI-JCE  ANALYSIS  «  12  O'CLOCK  POSITION 

.0225 


MODEL  3b  ENGINE 

1  . 

7.814  .05P 

0.0  0.5  1.0 

1.  .97  .13 

.0  .5 

0.627  0  ,t>9b 

TEST  «UN  t-OA  -  WE  1 
3208.  13000. 

18.05  1 « .  45 

12.09  351.6 

’3. 

TEST  HUH  61  A  ■  w*  T 
3266.  13000. 

17.45  1 o . 29 

15.16  351. 

33. 

TEST  HUN  62A  -  WET 
3295.  13000. 

14.82  1  «  ,'25 

12.11  400.2 

33. 

TEST  HUN  63m  -  wE f 
3259.  liOOO. 

17.4  1 1> .  30 

15.30  396. « 

TE^T  HIIn  64 A  -  wFT 
3222.  13000. 

17.40  16,35 

12.10  452.1 

TEST  HUN  65A  -  WK'I 
2831  .  1 3000 . 

18.49  16.21 

jj.35  452.1 

TEST  RUN  66 A  -  wET 
3253.  13000. 

1.02  20. 66 

J2.02  3  j  . . 4 

TEST  HUN  67A  -  wET 
3262.  13000. 

1,95  20.70 

15.17  351.5 

TEST  HUN  68A  -  WET 
3297.  13000. 

2:83  20,58 

12.03  399.5 

Te£t  HUN  69a  -  WET 
3280.  13000. 

3,53  20.65 

15.16  401.3 

tEST  RUN  70A  -  WET  CL'HPFCTED 
3266.  13000.  231.2 

2,01  20.7  .  .763  20.5 

12.10  451.1 

33, 

TEST  HUN  71 A  -  wET  CORRECTED 


0.125 

1.5  i.u  2.5 
1 .47-1 .30  -.94  -.81 
1.0  1.5  7.0 

0.445  0.897  n.609 

CORRECTED 
731.0 

.901  70.5 


C  P  H  H  (■.  C  T  F  D 
234.1 

.6^0  70.6 


CORRECTED  ' 
235,2 

.887  70.9 


cnx»F.riED 

234.5 

. o9o  20.6 


corrected 

233.0 

.694  20.7 


CUkHECI  El) 
234.7 

70.6 


CDhRc-CTFO 

231.0 

.762  70.5 


CORRECTED 

231.1 

.763  20.5 


corrected 

233.0 

.757  20.7 


CORRECTED 

232.2 

.760  20.6 


2.237 
3.0  3.5 


3291 

2,40 


13000. 

20,64 

452.4 


Mo 


.c 


•  •  V  , 

14.93 
33. 

TEST  RUN  72A  -  WEI  CORRECTED 
1637,  13000.  .  141.3 

-2,9?  23,76  1,944 

J|.06  350.6 


20.7 


5.4 


0.558  253. 

__  4.04.5  5.0 

-.77  -.81  -.92-1.09 
2.5  3.0  3.5 

0.317  0.328  0.350 


4.5 

0.350 
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.050 


8 • 5  Sample  Input  Data  (Cent.1! 
!Jm!T  “""hooo:  “f’.coHBtrrED 

jf:iT  IhVs 

4sfT  »F1  COKRKerPD 


S.5 


1645. 

-3.64 
11.90 
33, 

TEST  RIJN  75A 


13000. 

23.61 

401.3 


1  43 . 4 
1.917 


5.5 


1667. 

-4.66 
15.  J7 

fEST  RUN  76A 


13000. 

£3,54 

399.3 


1607. 

-2.61 
11.91 

Jest  run  77a 


13000. 

23,67 

452.1 


1552. 

T  4  .  r*  i 
15.40 

'Jest 

2389. 

;18.b 

12.14 
33. 


13000. 

?3.fc2 

452.8 


*fi  CORRECTED 

44.6 

l.*04  5.5 

WFT  CHkrecTFD 

142.6 

1.930  5. 5 

WF  i'  CORRECT  FD 
143.8 

1.917  5.5 


KUN1  C'ORRECTFD 


13000. 
24.21 
348. ft 


193.4 

1.442 


7.4 


?sfT  """SKS.-  *ET  ?SS»ECTKD 


;|lf{2  HTr  If  l,r 

lf-02  3bt;7  '-441 

-Jn!C  RUNi!°?s: 

?44j3  '-454 


:i9-43 

H-04 

4lfT  RUMJS„:  -ET  f™»§CTFD 


7.4 


6.5 


.  *  C. 

2541, 


7.5 


l5 

-J.TiT  *“1*11.* 

jj.K,  45IJ1  ‘-"J4  2.* 

l$V  "“Mooo: 

j|;?1  iU°.l  »&■» 
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APPENDIX  C  Table  C-l  Summary  of  Dry  Run  Skin  Temperature  Prediction  -  6  o'clock  Position 
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APPENDIX  C  Table  C-2  Summary  of  Dry  Run  Skin  Temperature  Prediction  -  9  o'clock  Posi_tjoni 
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